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ABSTRACT  OP  THE  DISSERTATION 


A  Perturbation  iheory  of  the  Heisenberg  Antiferromagnet 

by- 

Donald  Lynd  Bullock 

University  of  California,  Los  Angeles,  1963 
Assistant  Professor  David  I,  Paul,  Chairman 

The  determination  of  some  physical  parameters  of  the 
ground  state  for  the  Heisenberg  antiferromagnet  Is  con¬ 
sidered,  This  problem  has  been  completely  solved  hereto¬ 
fore  only  In  the  cases  of  the  energy  and  short  range  order 
parameters  for  the  linear  chain  with  spin  one-half.  In 
this  dissertation,  the  problem  of  more  general  lattices 
and  of  arbitrary  spin  is  considered,  and  the  long  range 
order  parameter  Is  treated  In  addition.  The  results  are 
reported  In  terms  of  series  expansions  generated  by  means 
of  a  modified  Rayleigh- Schro dinger  perturbation  theory, 
which  Is  proposed  and  developed  in  this  dissertation.  The 
modification  amounts  to  a  process  whereby  the  definition 
of  the  zero  order  Hamiltonian  is  changed  through  the  pre¬ 
diction  and  Inclusion  of  certain  infinite  classes  of  terras, 
whose  first  members  appear  In  the  original  perturbation 
series.  The  final  zero  order  Hamiltonian  obtained  is  the 
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Ising  model.  Explicit  expressions  for  the  ground  state 
parameters  are  given  through  fourth  order  for  the  linear 
chain,  plane  quadratic,  and  simple  cubic  lattices  with 
arbitrary  spin.  The  calculation  is  carried  through  to  de¬ 
termine  the  energy  series  through  six  orders  also  for  the 
linear  chain  with  spin  one-half,  A  comparison  with  experi¬ 
mental  determinations  of  the  long  range  order  in  real  anti- 
ferromagnets  is  made.  In  addition,  the  problem  of  the  first 
excited  states  is  considered.  A  perturbation  calculation 
similar  to  that  performed  for  the  ground  state  is  carried 
out  for  the  partition  function  of  the  Heisenberg  antiferro- 
magnet.  The  low  temperature  terms  of  this  expansion  are 
shown  to  be  related  to  the  energies  of  the  ground  state  and 
first  excited  states.  The  calculation  of  the  first  excited 
state  energies  is  carried  out  to  second  order.  A  complete 
analysis  of  the  relationship  of  this  perturbation  method 
for  the  ground  state  and  first  excited  states  with  the  spin 
wave  treatment  of  the  same  states  is  carried  out. 


I,  INTRODUCTION 


This  dissertation  concerns  some  physical  properties 
of  the  ground  state  of  the  Heisenberg  antif erromagnet.  The 
approach  is  one  of  developing  approximate  series  represen¬ 
tations  for  various  physical  parameters  by  means  of  a  mod¬ 
ified  Raylelgh-Schrodinger  perturbation  theory. 

For  classical  and  seml-classical  treatments  of  antl- 
ferromagnetlsm,  the  problem  of  the  ground  state  has  been 
largely  Ignored,  and  the  classical  antiferromagnetic  ground 
state,  or  Neel  state,  has  been  substituted  for  the  true 
ground  state.  Interest  In  the  true  ground  state  remains, 
however,  because  the  nature  of  the  problem.  Is  such  that 
its  solution  lies  in  the  direction  of  an  exact  solution  for 
the  entire  temperature  range.  Furthermore,  there  Is  the 
question  of  Just  how  much  similarity  there  Is  between  the 
physical  characteristics  of  the  true  ground  state  and  the 
often  substituted  Neel  state. 

In  Chapter  II,  the  Heisenberg  model  is  discussed,  with 
emphasis  on  the  problems  Inherent  In  treatments  of  the 
antiferromagnetic  lowest  energy  eigenstate.  Chapter  III 
Involves  the  derivation  of  an  effective  Hamiltonian,  which 
Is  subsequently  used  In  a  modified  perturbation  theory. 

This  modified  perturbation  theory  Is  then  proposed.  The 


3 


modification  amounts  to  a  process  whereby  the  definition 
of  the  zero  order  Hamiltonian  is  changed  through  the  pre¬ 
diction  and  Inclusion  of  certain  Infinite  classes  of  terms 
whose  first  members  appear  In  the  original  perturbation 
series.  The  final  zero  order  Hamiltonian  obtained  In 
Chapter  III  Is  the  Ising  model.  Perturbation  series  are 
then  reported  through  the  fourth  order  for  the  energy^ 
short  range,  and  long  range  order  parameters  for  typical 
lattices  of  one,  two,  and  three  dimensions  and  arbitrary 
spin.  In  addition,  an  analysis  of  the  relationship  of 
these  perturbation  corrections  to  the  corrections  gener¬ 
ated  by  the  spin  wave  theory  Is  carried  out.  In  Chapter 
IV,  the  same  treatment  Is  applied  to  the  first  excited 
states  of  the  Heisenberg  antiferromagnet,  along  with  a 
subsequent  analysis  of  Its  relationship  to  the  spin  wave 
treatment  of  the  same  states. 

The  Importance  of  this  particular  quantum  mechanical 
treatment  Is  that  it  is  a  new  treatment  of  the  Heisenberg 
antiferromagnet  at  low  temperatures.  Further,  this  new 
method  Is  analyzed  and  then  compared  with  the  standard 
quantimi  mechanical  approximation  technique,  l.e,,  the 
spin  wave  theory,  thus  giving  additional  insight  Into  the 
problem. 


II.  THE  HEISENBERG  MODEL  AND  THE 
GROUND  STATE  PROBLEM 

A,  Discussion  of  the  Model 

Neel^  was  the  first  to  propose^  and  extensively  Inves¬ 
tigate,  the  antiferromagnetic  state.  The  basis  for  his 

treatment  of  this  state  Is  the  model  Hamiltonian  introduced 
2 

by  Heisenberg  in  his  first  paper  on  the  theory  of  ferro¬ 
magnetism.  This  model  continues  to  be  the  most  common 
starting  point  for  treatments  of  cooperative  magnetic 
phenomena  In  both  ferro-  and  antiferromagnetism,  as  well 
as  in  the  intermediate  case  of  ferrlmagnetlsm.^ 

The  Heisenberg  model  is  useful  In  that  It  seems  to 
be  sufficiently  representative  of  the  physics  of  coopera¬ 
tive  magnetism  to  predict  with  some  consistency^  the  re¬ 
sults  of  experiment.  In  spite  of  this  success,  however, 
there  remain  Important  theoretical  problems  with  regard 
to  the  derivation  and  Interpretation  of  the  model.  There 
Is  always  the  need  to  derive  the  model  unambiguously  from 
the  full  Hamiltonian  expressed  In  terras  of  electron  and 

nuclear  coordinates,  or  to  show  that  something  quite  like 

4 

It  exists.  Such  a  derivation  would  show  clearly  what 
constitute  the  limitations  and  omissions  of  the  model.-' 
There  Is  also  the  need  for  the  complete  mathematical 
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solution  of  the  model  Itself ^  i.e„j(  determination  of  the 
eigenvalue  spectrum  and  degeneracy  scheme.  The  emphasis 
in  this  dissertation  is  with  regard  to  a  problem  which 
lies  In  the  latter  area. 

Complete  descriptions  of  the  eigenvalue,  spectrum  and 
degeneracy  scheme  for  the  Heisenberg  model  do  not  exists 
except  for  the  treatment  of  the  spin  one -half  coupled 
linear  chain  by  Bethe.  Even  for  that  case,  determina¬ 
tion  of  the  eigenvalue  spectrum  and  degeneracy  Is  given 
only  In  principle.  The  set  of  coupled  transcendental 
equations  which  result  has  been  solved  for  the  energy  of 
the  antlferroraagnetlc  ground  state”^  and  of  the  antiferro- 

O 

magnetic  first  excited  states,  as  well  as  for  some  low- 
lying  ferromagnetic  states.  The  great  bulk  of  states 
which  contribute  to  the  thermodynamical  properties  at 
finite  temperatures  remains,  however,  uninvestigated, 

Q 

Progress  with  exact  solutions  of  even  the  Ising  model, ^ 
which  Is  a  considerably  simplified  abbreviation  of  the 
Heisenberg  model,  has  so  far  failed  of  any  exact  three 
dimensional  treatment.  The  present  effort  is  directed 
toward  the  problem  of  the  antiferromagnetic  ground  state 
for  general  spin  and  typical  one,  two,  and  three  dimen¬ 
sional  lattices  from  a  perturbation  theoretic  point  of 
view. 

The  mathematical  statement  of  the  Heisenberg  model 
Incorporating  the  nearest  neighbor  approximation  Is, 
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H  =  -2J  I  (1) 

<JlO  ^  ^ 

where  J  Is  the  exchange  Integral,  and  the  sum  is  carried 
out  over  the  scalar  product  of  all  nearest  neighbor  pairs 
of  spin  angular  momentum  operators  associated  with  a  given 
lattice  structure .  Eq.  (l)  has  both  a  classical  and  a 
quantum  mechanical  Interpretation.  The  determinations  of 
the  lowest  energy  state  when  Sj  and  Sj^  are  considered  to 
be  ordinary  vectors,  regardless  of  whether  J  Is  less  than 
(antiferromagnetism)  or  greater  than  zero  (ferromagnetism), 
is  a  trivial  problem,  Por  J  >  0,  all  scalar  products  must 
be  maximized,  giving  rise  to  a  ground  state  configuration 
of  perfect  alignment  which  is  fully  determined  except  for 
a  directional  degeneracy^*^  arising  from  the  complete  Iso¬ 
tropy  of  (l).  Ordinarily  this  directional  degeneracy  is 
resolved  by  Introducing  a  finite  anisotropy  Into  Eq.  (l). 
The  directional  degeneracy  is  removed  completely  for  ferro¬ 
magnetism  by  the  Introduction  of  an  external  magnetic 
field.  For  antiferromagnetism  it  is  reduced  to  a  double 
degeneracy  by  the  Internal  crystalline  anisotropy. 

In  the  antiferromagnetic  case,  J  <  0,  the  minimization 
of  H  results  In  the  minimization  of  the  Individual  scalar 
products  only  if  the  total  lattice  can  be  divided  Into  two 
equivalent  Interlocking  sublattices  such  that  all  nearest 
neighbors  of  a  given  sublattice  point  lie  on  the  other 
sublattice. .  Again,  this  Is  fully  determined  except  for 


a  directional  degeneracy.  (In  Figure  1  are  shown  the 
classical  ground  state  configurations  for  typical  anti- 
f erromagnets , ) 

Quantum  mechanics  presents  no  special  problem  In  the 
determination  of  the  ground  state  for  ferromagnetism. 

It  Is  simply  the  quantum  mechanical  analog  of  the  classi¬ 
cal  ground  state^  i.e.j,  individual  lattice  point  spin 
systems  In  states  of  maxlmiun  angular  momentiun  along  some 
common  axis.  For  antiferromagnetism^  however^  the  deter¬ 
mination  of  the  quantum  mechanical  ground  state  is  a  prob¬ 
lem  of  considerable  difficulty.  In  faet^  only  certain 
parameters  of  the  antiferromagnetic  ground  state  have  been 
investigated. 

For  classical  and  seml-classical  treatments  of  anti- 
ferroraagnetlsm,  the  problem  of  the  groxmd  state  has  been 
largely  ignored,  and  the  classical  antiferromagnetic  ground 
state,  or  Neel  state,  has  been  substituted  for  the  true 
ground  state.  Interest  in  the  true  ground  state  remains, 
however,  because  the  nature  of  the  problem  Is  such  that  Its 
solution  lies  In  the  direction  of  an  exact  solution  for  the 
entire  model,  Techniques  which  successfully  treat  the 
ground  state  may  be  applicable  to  certain  other  states,  and 
eventually  inay  be  applicable  to  all  other  states.  Further¬ 
more,  there  Is  the  question  of  Just  how  much  similarity 
there  Is  between  the  physical  characteristics  of  the  true 
ground  state  and  the  often  substituted  Neel  state. 


B,  Theoretical  Treatments  of  the  Ground  State 


The  one  exact  treatment  of  an  antiferromagnetic  ground 
state  has  been  previously  mentioned.^  This  treatment  of 
the  spin  one -half  coupled  linear  chain  provides  a  criterion 
for  the  comparison  of  all  approximate  treatments.  Approxi¬ 
mate  methods  which  treat  well  the  ground  state  of  the  line¬ 
ar  chain  are  likely  to  treat  well  the  ground  state  for 
other  spins  and  lattices  If  the  method  Is  applicable. 

During  the  past  ten  years  there  has  been  much  activi¬ 
ty  related  to  the  determination  of  the  physical  character¬ 
istics  of  the  true  ground  state  by  means  of  approximation. 
The  first  quanttun  mechanical  treatment  for  arbitrary  spin 

systems  In  an  Interlocking  sublattice  structure  was  given 
11 

by  Anderson.  In  it  he  applied  what  was  qualitatively  a 
spin  wave  approximation  to  Eq.  (1).  The  results  are  valid 
for  large  values  of  the  total  spin  of  Individual  lattice 
point  systems  (all  lattice  point  systems  having  a  common 
total  spin)  and,  when  extrapolated  to  values  of  the  spin 
of  order  unity,  give  physically  Interesting  predictions  of 
the  sublattice  magnetization  (long  range  order).  Of 
course,  the  ground  state  eigenvalues  were  also  determined, 
but  since  they  are  not  susceptible  to  experimental  verifi¬ 
cation,  their  importance  Is  relegated  to  the  role  of  com¬ 
parison  with  the  eigenvalues  determined  by  other  methods. 
Anderson  did  point  out,  however,  that  In  ail  cases  treated 
his  approximate  eigenvalue  lay  within  rigorously  determined 


12 

limits  for  the  true  eigenvalue — thus  establishing  that 
the  method^  at  the  very  leasts,  treats  the  energy  adequate¬ 
ly. 

A  special  treatment  of  the  short  range  order  for  the 

spin  one -half  coupled  linear  chain  was  considered  by 

Hulthen.'^  The  variational  technique  which  he  used  was 

1^ 

improved  upon  by  Kasteleijn  and  extended  to  an  aniso¬ 
tropic  Hamiltonian  which  is  a  function  of  a  parameter 
and  has  as  its  limiting  cases  the  Heisenberg  and  Ising 
models,  Kasteleijn  also  calculated  the  long  range  order 
for  this  case,  and  the  results  obtained  in  the  Heisenberg 

limit  were  qualitatively  the  same  as  the  predictions  of 

l4 

the  spin  wave  theory  of  Anderson,  However;,  Orbach  was 
able  to  calculate  the  energy  and  short  range  order  for 
this  case  as  a  function  of  the  anisotropy  without  approxi- 
matioH;  and  his  results  are  materially  different  from 
those  of  Kasteleljn,  This  cast  serious  doubt  upon  the 
reliability  of  the  variational  treatments  as  applied  to 
the  determination  of  antiferromagnetic  ground  state  param¬ 
eters. 

Prior  to  Orbach 's  discovery,  the  variational  treat¬ 
ment  of  Kasteleljn  had  been  extended  to  other  cases, 

Taketa  and  Nakaraura^^  extended  the  Kasteleljn  treatment, 
including  the  anisotropy  parameter,  to  two  and  three  di¬ 
mensional  lattices.  Their  results  in  the  isotropic  limit 
were  in  disagreement  with  those  of  Anderson  for  these 

g 


cases ^  since  their  long  range  order  parameters  exhibited 
a  behavior  as  a  function  of  the  anisotropy  which  was  simi¬ 
lar  to  that  of  Kasteleijn's  for  the  linear  chaino 

Marshall, using  the  isotropic  (Heisenberg)  Hamil¬ 
tonian  and  a  combination  of  the  variational  method  and  the 
Bethe-Pelrls  approximation,  ‘  calculated  the  energy,  short 
range,  and  long  range  order  parameters  for  two  and  three 
dimensional  spin  one-half  systems,  obtaining  predictions 
which  were  also  in  conflict  with  the  spin  wave  theory.  On 
the  other  hand,  Kubo*^  used  a  combination  of  the  varia¬ 
tional  method  and  the  spin  wave  approximation  itself  to 
treat  lattices  of  one,  two,  and  three  dimensions  and  ar¬ 
bitrary  spin,  obtaining  results  which  were  in  agreement 
with  the  spin  wave  theory  except  in  the  case  of  one  dimen¬ 
sion,  For  the  one  dimensional  lattice,  Kubo  predicted  an 
ordered  ground  state,  a  result  which  is  in  contradiction 
to  the  spin  wave  theory. 

C  o  Recent  Perturbation  Treatments  of  the 
Ground  State 

In  the  past  few  years,  three  new  perturbation  treat¬ 
ments  of  Eq,  (l)  for  the  ground  state  have  been  introduced. 
The  original  perturbation  treatment  of  Eq,  (l)  consisted 
of  a  division  into  the  spin  wave  Hamiltonian  plus  a  per¬ 
turbation,  Walker^^  reports  a  perturbation  treatment  of 
the  energy,  short  range  order,  and  long  range  order  of  the 
spin  one -half  coupled  linear  chain  through  sixth  order. 
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using  the  Islng  model  as  a  zero  order  Hamiltoniano  For 
comparison  with  the  exact  treatments  of  the  energy  and 
short  range  order,  he  derives  a  series  solution  expressed 
in  terms  of  the  same  parameter  as  was  used  in  the  pertur¬ 
bation  treatment  dl.rectly  from  the  integral  equations  as 
given  by  Orbach,  The  agreement  between  the  two  series 
is  complete  through  sixth  order. 

The  division  of  Eq„  (l)  used  by  Davis  to  develop  a 
"linked  cluster"  expansion  for  the  energy,  short  range 
order,  and  long  range  order  is 

Hq  =  -2J^I^[S^SJ-(S-S^)(S+S^)]  (2) 

=  “2J^£^[(S-S^)(S+s|)  +  (l-a)(sjsg+sysg)]. 

Although  does  not  retain  as  much  of  Eq„  (l)  as  does  the 
Ising  model,  it  describes  a  system  of  "independent  parti¬ 
cles,  "  for  which  Davis  is  able  to  prove  the  existence  of 
a  "linked  cluster"  expansion.  Using  this  technique,  he 
is  able  to  report  perturbation  corrections  through  the 
seventh  order  in  for  all  three  ground  state  parameters. 
In  order  to  improve  convergence,  he  also  reports  contribu¬ 
tions  from  higher  orders,  which  he  is  able  to  predict  on 
the  basis  of  regularities  evident  in  the  first  seven  opders. 
He  obtains  energy  eigenvalues  comparable  with  those 


obtained  by  other  workers  and  predicts  order  In  the  ground 
state  for  the  Heisenberg  limit  of  Eq.  (2)  for  all  three 
dimensions o 
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The  zero  order  Hamiltonian  considered  by  Roth  is 
the  long  wavelength  limit  of  a  Fourier  transform  of  Eq„ 

(1)5  Isfioj 


Ho  =  (3) 


Hn  =  -2jI'^^[s5  ]. 


The  prime  on  the  summation  symbol  indicates  the  omission 
of  the  long  wavelength  term^  K  =  0.  The  transformed  oper¬ 
ators  are  of  the  form 


X  ?  V2r 
Sj  =  (f)  I  exp  (IK.J)S- 
^1  ^  J 


,x 

r 


(4) 


(-lK.k)Sj^. 


The  sum  over  J  represents  a  sum  over  lattice  sites  be¬ 
longing  to  one  of  the  sublattlceSj  and  the  sum  over  k  rep¬ 
resents  a  sum  over  sites  belonging  to  the  other  sublattice. 
N  is  the  total  number  of  lattice  sites  in  the  crystal,  and 
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Y„  =  Z  exp  (iK->(5)  (5) 

it  Q 

with  5  representing  the  vector  separation  between  nearest 

neighbors.  The  advantage  of  this  particular  division  of 

Eq„  (1)  Into  zero  order  and  perturbation  parts  is  that  Hq 

is  Isotropic.  The  convergence  properties  of  the  series 
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are  similar  to  the  Davis  series  for  the  energy  before  the 
inclusion  of  the  terms  from  orders  higher  than  the  seventh. 
No  short  range  order  or  long  range  order  series  have  been 
reported  using  Eq„  (3), 

In  addition  to  the  above  three  perturbation  treatments 
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of  the  ground  state.  Boon  has  developed  an  energy  series 
based  on  the  Islng  model  as  zero  order  Hamiltonian  which 
Is  arbitrary  In  the  spin  and  dimension  of  the  lattice. 

His  results,  however,  are  different  from  the  results  ob¬ 
tained  later  In  this  dissertation  for  the  general  case. 

They  differ  with  Walker  ^  also  for  the  linear  chain. 

Table  I  gives  the  qualitative  results  with  regard  to 
the  prediction  of  a  finite  sublattice  magnetization  for 
all  of  the  treatments  mentioned  previously  which  report 
such  a  quantity.  The  results  are  seen  to  be  quite  con¬ 
tradictory,  but  In  general,  the  variational  treatments 
predict  zero  sublattice  magnetization  and  the  perturba¬ 
tion  treatments  predict  a  finite  sublattice  magnetization. 
The  explanations  generally  given  for  this  situation  are 
the  following;  (1)  the  perturbation  treatments  begin  with 


the  perfectly  ordered  Heel  state  and  then  carry  the  per¬ 
turbation  through  only  a  few  orders  in  the  calculation^, 
resxiltlng  in  a  Mas  towards  predicting  sublattice  mag¬ 
netization  in  the  ground  state^,  and  (2)  the  variational 
treatments  reflect  the  existence  of  a  statistically  large 
number  of  low-lying  states  which  exhibit  zero  sublattice 
magnetization.  Thus  the  existence  of  a  sublattice  mag¬ 
netization  in  the  antiferromagnetic  ground  state  seems  to 
be  unanswered  in  spite  of  the  numerous  attempts  to  deter¬ 
mine  it. 

From  an  experimental  point  of  vlew^,  real  antiferro- 

magnets  do  exhibit  a  sublattice  magnetization  at  low  tem- 
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peratures,  If  the  Heisenberg  model  is  to  give  a  cor¬ 

rect  picture  for  antiferromagnetism,  then  it  must  predict 
a  finite  sublattice  magnetization  for,  at  the  very  least, 
the  ground  state  region.  If  the  Heisenberg  model  is  as¬ 
sumed  to  be  essentially  correct  for  antiferromagnetism, 
then  it  must  be  concluded  that  the  number  of  low-lying 
states  exhibiting  zero  sublattice  magnetization  do  not 
statistically  outnumber  those  exhibiting  a  finite  sub¬ 
lattice  magnetization.  This  being  the  case,  the  predic¬ 
tions  of  the  variational  treatments  cannot  be  explained 
away  in  manner  previously  suggested.  Of  course,  if  the 
Heisenberg  model  does  not  give  the  picture  for  a  real 
antlferromagnet,  then  it  may  be  true  that  the  low-lying 
states  exhibiting,  zero  sublattices  magnetization  outnum- 


ber  those  with  finite  sublattice  magnetization  for  the 
model o 

Recently there  have  been  some  new  speculations  on  the 
existence  of  a  sublattice  magnetization  for  the  Heisenberg 
model  antiferromagnetic  ground  state  which  affect  the  in¬ 
terpretation  of  both  the  perturbation  and  variational  ap- 
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proaches,  Pratt  has  pointed  out  that^  if  the  ground 
state  is  non-degenerate^  the  time  reversal  symmetry  of  the 
Hamiltonian  of  Eq.  (l)  requires  the  sublattice  magnetiza¬ 
tion  to  be  zero.  In  addition^  the  same  arguments  may  be 
applied  to  anisotropic  Hamiltonians  of  the  form  of  Eq,  (2), 
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KaraylaniSji  et  ad.;,  have  shown  that  the  ground  state  of 
the  spin  one-half  coupled  linear  chain  is  non-degenerate 
for  arbitrary  value  of  the  anisotropy  (except  for  the  Islng 
limit) j  so  that  the  arguments  of  Pratt  certainly  apply  to 
this  case.  In  the  Ising  limits  there  is  a  degeneracy  in 
the  ground  state  consisting  of  the  two  time  reversal  sym¬ 
metric  Neel  states  for  the  linear  chain  which  are  shown  in 
Figure  1,  Each  of  these  states  alone  exhibits  perfect  sub- 
lattice  magnetization^,  but  the  correct  linear  combination 
of  them  which  fits  onto  the  ground  state  wave  function  near 
the  Islng  limit  is  the  normalized  siun  of  the  two  Neel 
states,  which  exhibits  zero  sublattice  magnetization.  This 
makes  the  prediction  of  zero  sublattice  magnetization  for 
the  anlso tropically  coupled  linear  chain  with  spin  one-half 
consistent  throughout  the  range  of  the  anisotropy  parameter. 


In  view  of  this  development  for  the  ground  state  of 
the  spin  one -half  coupled  linear  chaln^,  It  is  Important  to 
consider  what  effect  a  similar  development  for  other  cases 
may  have  on  the  interpretation  of  sublattice  magnetization 
calculations,  such  as  the  one  performed  In  this  disserta¬ 
tion.  Since  a  sublattice  magnetization  must  be  accepted 
as  a  property  of  all  real  antlferromagnets  at  very  low 
temperatures,  and  since  the  Heisenberg  model  seems  essen¬ 
tially  correct  for  the  discussion  of  many  properties  of 
real  antlferromagnets,  it  will  be  assumed  without  further 
discussion  that  the  time  reversal  symmetry  argument  is 
disqualified  in  some  way  for  the  purposes  of  this  disser¬ 
tation,  This  might  come  about  through  the  degeneracy  of 
the  ground  state  for  the  important  three  dimensional  cases, 
or  through  the  destruction  of  the  time  reversal  symmetry 
by  some  kind  of  non- time-reversal -symmetric  interaction. 
Thus  a  finite  sublattice  magnetization  will  be  assumed 
possible,  and  the  Neel  state  will  be  chosen  for  a  zero 
order  wave  function. 

We  now  turn  to  the  main  purpose  of  this  dissertation, 
which  is  to  propose  a  new  perturbation  approach  to  the 
physical  properties  of  the  antiferromagnetic  ground  state. 


Ill,  PERTURBATION  THEORY  OP  THE  ANTIPKRROMGNETIC 

GROUND  STATE 

Ao  The  Effective  Hamiltonian 

1.  Division  of  the  Hamiltonian 

In  this  chaper  we  propose  a  new  perturbation  ap¬ 
proach  to  the  problem  of  the  antiferromagnetic  ground 
state.  We  generate  a  series  for  the  energy  which  Is  gen¬ 
eral  In  the  spin  and  valid  for  any  interlocking  sublattice 
system.  The  Islng  model  is  taken  to  be  the  zero  order  ap¬ 
proximation,  and  In  this  way  our  series  differs  from  both 
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the  spin  wave  theory,  the  linked  cluster  perturbation 
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treatment,  and  the  diagram  technique.  Such  a  zero  or¬ 
der  approximation  affords  a  means  of  determining  the 
ground  state  parameters  without  being  subject  to  such  crit¬ 
icisms  as  the  Inclusion  of  non-physical  states  and  the 
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omission  of  the  dynamical  Interaction,  Also,  the  uncon¬ 
ventional  way  in  which  we  obtain  the  Islng  model  for  zero 
order  approximation  affords  a  means  of  determining  the  sub¬ 
lattice  magnetization  directly  from  the  energy  series  by 

29 

means  of  the  Feynman  theorem.  ^  This  is  not  a  property  of 
an  energy  series  based  strictly  on  the  Islng  zero  order 
approximation , 

The  method  outlined  below  is  most  similar  to  the  work 
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Pn 

of  Davis j  but  it  di.ffers  in  two  important  respects. 

Firsts  it  does  not  require  the  use  of  a  linked  cluster 
criterion  for  the  determination  of  physical  contributions 
to  the  perturbation  series.  Second,  it  leads  to  the  com¬ 
plete  recovery  of  the  Islng  model  as  zero-order  Hamiltonian 
through  four  orders  of  the  perturbation  for  the  general 
case,  and  six  orders  for  the  spin  one-half  coupled  linear 
chain.  This  second  difference  is  probably  the  more  sig¬ 
nificant  since  it  affords  a  basis  upon  which  an  instruc¬ 
tive  comparison  with  the  spin  wave  treatment  of  the  ground 
state  is  possible, 

2 ,  Discussion  of  the  Anisotropy 

We  begin  our  perturbation  treatment  vflth  the  con¬ 
struction  of  an  effective  Hamiltonian  to  replace  Eq.  (l)„ 

We  first  write  Eq,  (1)  in  a  slightly  different  form,  l,e, 

H  =  2J  I  S^SJ  +  J(l-a)  I  (S^sf+SM) ,  (6) 

Introducing  the  parameter  a  and  allowing  this 
parameter  to  take  on  non-zero  values,  we  are  modifying 


* 

Here  we  have  written  the  Hamiltonian  so  that  J  is 
positive  definite  for  antiferromagnetism,  and  thus  J 
carries  no  Implicit  sign  change  for  any  expression  involv¬ 
ing  J, 


Eq„  (l)  to  include  the  effects  of  anisotropy „  In  a  physi¬ 
cal  sense^  we  are  introducing  a  form  of  the  anisotropic 
exchange  interaction  which  was  derived  by  Van  VLeck'^  for 
exchange  coupled  spin  one-half  magnetic  ions.  This  inter¬ 
action  takes  the  following  form  for  arbitrary  spin. 


H. 


Int. 


J  1^' 


(7) 


where  J  and  k  refer  to  different  lattice  sites ^  and^  and 
V  denote  the  three  components  of  spin  angular  momentumo 
The  may  be  determined  from  a  second  order  perturba¬ 

tion  calculation  involving  the  spin-orbit  coupling  and 
wave  function  overlap  between  nearest  neighbor  lons^  For 
S  =  1/S!j  all  with  /a  ^  V  are  identically  zero.,  giving 

rise  to  a  form  which  is  comparable  with  the  anisotropy  in¬ 
troduced  in  Eq,  (6),  Even  for  cases  in  which  S  1/2^  our 
form  of  the  anisotropy  is  comparable  with  Eq„  (7)  provided 
the  for are  sufficiently  smalls 

In  a  mathematical  sensej,  the  introduction  of  the 
parameter  a  allows  Eq.  (6)  to  represent  either  the  Islng 
models  (a=l)j,  or  the  Heisenberg  models  (a=0)„  In  addition^, 
for  a  perturbation  calculation  based  upon  the  Islng  model 
as  zero  order  approximation^  the  factor  (1-a)  serves  as  a 
perturbation  parameter. 
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3.  Transformation  to  Reciprocal  Lattice  Space 

The  commutation  relations  for  spin  operators  at¬ 
tached  to  the  various  lattice  points  are 


-  V*y 


(8) 


and 


Where  we  have  Introduced  the  usual  relationship  between 
spin  components  and  the  raising  and  lowering  operators^ 


S'*'  =  +  IS^, 


(9) 


s"  -  -  IS^. 

Using  these  relationships  and  the  transformation  to  recip¬ 
rocal  lattice  space,  l.e.. 


^  -  (2/N)^/^  I  exp  (lA*J)S. 


(10) 


=  (2/N)^/^  I  exp  (-lA*k)S^ 


and 
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Where  the  summations  over  J  and  k  represent  summations  over 
the  lattice  vectors  of  the  two  raspectlve  sublattices,  we 
find  for  Eq.  (6), 

where 

■2'x  =  5  (l/\‘5).  (12) 

6 

The  vector  S  represents  the  separation  of  nearest  neighbor 
lattice  sites.  The  conanutatlon  relations  In  reciprocal 
lattice  space  are 

(13) 

and 

We  define  our  Neel  state  (l.0>)  to  be 

st  1  0>  =  0, 


(14) 


^jo>  -  (S/S)^''^Si^lO>, 


%l  0>  .  0, 

and 

|  0>  =  -  (N/2)^/^S5;^  |o>. 
2  0 


Thus  10>  corresponds  to  a  product  of  spxn  eigenfunctions^, 
one  for  each  lattice  site^  in  which  all  eigenfunctions  on 
the  J  sublattice  (henceforth  referred  to  as  sublattice  #1) 
have  maximum  z  component  of  spin,  and  all  those  on  the  k 
sublattice  (henceforth  referred  to  as  sublattice  #2)  have 
minimum  z  component  of  spin.  This  is  perfect  (l.e.,  clas¬ 
sical)  anciferromagnetic  ordering. 

Next,  we  express  those  states  in  which  the  z  com¬ 
ponent  of  spin  has  been  reduced  by  one  on  sublattice  #1, 
and  Increased  by  one  on  sublattice  #2,  We  do  not  need 
to  investigate  states  in  which  the  total  z  component  of 
spin  for  the  whole  lattice  is  non-zero,  since  the  Neel 
state  has  zero  z  component  of  spin  and  Eq,  (11)  conserves 
the  total  z  component  of  spin.  The  states  described  above 
are  of  the  form 


where  J  and  k  are  arbiti'ary  lattice  sites  in  each  of  the 
two  sublattices.  In  reciprocal  lattice  space  Eq,  (15)  be¬ 
comes 

(2/N)  I  exp  [-l(A*J+M*k)3sT  st  1.0>„  (16) 

The  natural  extension  of  Eq.  (15)  to  other  states  with 
total  z  component  of  spin  equal  to  zero^,  e,g,,  for  a  typi¬ 
cal  state 


3;  )]  X  [(sj  )(sj 

^n  ^1  ^ 


)]lG>, 

'  (17) 


results  in  a  reciprocal  lattice  wave  function  which  is  a 
linear  combination  of  terms  of  the  form 


where  Za^  «=  take  Iji  as  represen¬ 

tative  of  a  typical  wave  function  involving  the  reciprocal 
lattice  spin  operators  and  having  total  z  component  of 
spin  equal  to  zero. 


Construction  and  Physical  Sij: 
Effective  Hamiltonian 


liflcance  of  the 


We  now  define  an  effective  Hamiltonian  to  replace  Eq 
(11),  We  begin  by  writing  Eq,  (11)  as  the  sum  of  three 


terms 


«z  -  (19) 

A  12 

A  12 

■Me  '  J(i-^)I*aMa  MA  . 

A  1  c 

The  second  and  third  terms  of  Eq.  (19)  have  the  following 
physical  significance;  the  second  term  decreases  the  z 
component  of  angular  momentum  on  sublattice  #1  by  one  unit 
while  simultaneously  raising  the  z  component  of  angular 
momentum  on  sublattice  #2  by  one  unit.  When  applied  to 
the  Neel  state  described  by  Eq.  (14),  we  may  think  of 
as  raising  the  Neel  state  to  an  excited  state.  This  prop¬ 
erty  of  Hy  Is  the  reason  for  choosing  the  subscript  U  to 
denote  this  term  of  the  Hamiltonian,  The  operator  H-g 
raises  a  given  state  "up"  from  the  point  of  view  of  the 
n\imber  of  excitations  present  over  the  original  Neel  state 
Correspondingly,  the  operator  Hg  lowers  the  number  of  ex¬ 
citations  "down"  from  the  number  of  excitations  present. 

We  shall  leave  Hg  and  Hp  unchanged  from  their  form  In  Eq, 
(19),  and  shall  Incorporate  them  directly  into  our  new 
effective  Hamiltonian.  We  derive  this  form  by  applying 

to  the  typical  wave  function  of  Eq,  (l8).  H  may  then  be 

z 
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commuted  through  all  of  the  operators  in  until  it  oper¬ 
ates  directly  on  the  Neel  state.  The  resulting  extra  terms 
all  Involve  z  component  spin  operators  which  are  also  com¬ 
muted  through  all  operators  until  they  operate  upon  the 
Neel  state.  Since  the  Neel  state  is  an  eigenftinction  of 
all  z  component  spin  operators^  we  may  replace  them  with 
their  eigenvalues  to  obtain 


1. 

Eq.  (20)  indicates  that  the  Islng  Interaction  may  be 
expressed  by  the  sum  of  two  parts.  The  first  part  consists 
of  a  Hamiltonian  which  describes  Independent  excitations 
above  the  Neel  state,  and  the  second  part  describes  an  in¬ 
teraction  between  these  excitations.  If  we  write  the 
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Hamiltonian  of  the  interaction  as  then  our  effective 

■Rl 

Hamiltonian  takes  the  form*^ 


The  operator  fom  for  Hq  is 


but  we  may  always  represent  Eq.  (22)  by  its  eigenvalues 
in  the  work  that  follows. 

Having  defined  the  effective  Hamiltonian,  we  are  now 
in  a  position  to  apply  perturbation  theory  to  Eq.  (21). 


B.  The  Application  of  a  Modified  Rayleigh-Schro dinger 
Perturbation  Theory  to  the  Neel  State 

1.  The  Modified  Perturbation  Method 

We  divide  Eq.  (21)  into  zero  order  and  perturba¬ 
tion  Hamiltonians, 


Hq  “  Eq  2^  ^ 

A  ' 


+  Hj  +  Hp 


The  (|J  as  defined  constitute  eigenfunctions  of  Hq. 
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We  do  not  consider  a  particular  choice  of  to  span  the 
(2S+1)^  dimensional  spin  space  of  Interest,  but  only  refer 

•31 

to  the  fact  that  such  a  choice  exists.-^  Taking  the  Neel 
state  for  our  zero  order  wave  function,  l|3  q,  we  have  upon 
application  to  the  Raylelgh-Schrodlnger  perturbation 
theory, 

•*0^0  =  ^^0'  (25a) 

+  Hjf  0  “  1  =liPo>  (25’=) 

«q92  “1^1  -  '’of  2  +  1  +  '’sf  0'  (25°) 

etc. 

The  first  order  correction  to  the  wave  function,  (p 
may  be  expressed  as  a  sum  of  contributions  Involving  dif¬ 
ferent  nxunbers  of  "independent"  excitations,  i.e., 

f  1  ■  ^  •••'  (25) 

= 

n 

where 

Tl^^  =  aj°Ho>, 
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etCo  Here  we  have  restricted  to  contributions  of  the 
form  of  Eqo  (l8)  having  an  equal  number  of  spin  operators 
for  each  sublattice.  Such  contributions  satisfy  the  eigen¬ 
value  equation. 


'total  T1 


where  operator  for  the  total  z  component  of 

spin  for  the  entire  lattice,  l.e.,  ^total  ~  ^  ^ k 
That  only  such  contributions  are  necessary  for  the  expan¬ 
sion  of  (or  any  higher  order  correction  to  the  wave 
function  )  follows  from  the  two  properties  of  and  Cp  q, 

^totalH^O  “  ^^total'’  %)  “  0.  In  other  words,  the 

perturbation  H-j_  does  not  connect  (pQ  with  any  states  having 
a  different  total  z  component  of  spin.  We  have  previously 
discussed  this  property  with  respect  to  the  total  Hamil¬ 
tonian,  Eq.  (ll),  but  it  is  necessary  here  to  specify  it 
for  the  individual  zero  order  and  perturbation  Hamiltonl- 


Thus  we  have  expressed  in  terms  of  a  set 
where  is  a  linear  combination  of  states  of  the  lat¬ 

tice  wherein  there  are  n  spin  operators  Si  on  sublattice 

+ 

#1,  and  n  spin  operators  S  ^  on  sublattice  #2,  We  notice 
(  )  '^2 

that  the  pp  ■'  are  orthogonal  in  n.  Further, 


=  (Eg  +  nA)p^''), 


(27) 


and  thus  from  Eq.  (25b)  we  get 


a  Icp(“)  .  (Ej^-Hj,)cpg,  (28) 

Since  =  0^  and  H^^CpQ  =0.  On  the  left  of  Eq.  (28) 

are  a  set  of  contributions  to  arising  from  orthogonal 
subspaces  of  the  original  (2S+l)N  dimensional  spin  space, 
and  Eq,  (28)  must  be  satisfied  separately  In  each  of  these 
subspaces.  The  components  of  have  three  Important 
properties  with  respect  to 

(29) 

and 

Hy  transforms  a  function  defined  entirely  In  the 

subspace  n.  Into  a  function  which  Is  defined  en¬ 

tirely  In  the  subspace  n  +  1,  Hj  leaves  unchanged  the  sub¬ 
space  in  which  Is  defined,  and  transforms  cp 

Into  a  function  0^^  which  Is  entirely  defined  within 
the  sub space  n  =  1. 

In  the  n  =  0  subspace  we  have  from  Eq,  (28), 
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and  therefore 


0 


^  0  O 


For  n  =  1  we  have 


-  Hgipo, 

and  therefore 

-  -  i  •  %to- 

Por  n  >  2  we  find  that  0,  so  the  first  order  of  the 

perturbation  yields 


E  “  V  (30) 

and 

since  =0  In  general. 

Extending  this  procedure  to  higher  order and  intro¬ 
ducing  the  simplified  notation 


(31) 


we  find  for  the  energy  corrections  through  the  sixth  order 


\  «  0,  (32) 

Eo  =  -  ^KDU><, 

2  A 

E3  =  1/a2<I>IB>, 

E4  a  1/a3[<DU>^  -  <DI^U>  -  ~<D^U%], 

E3  =  1A^[-3<DBX.BITX>  +  <DI^U>  +  <D^UIU>  + 

Eg  =  l//\^[-2<m>^  +  4<DUXDIU>  +7/4<DUXD^uS 

+  2<DH}>^  "  ^IDIU^>  -  <D^UI^U> 

122  4  1222 

-  ^D'^UDU  >  -  <2>I^IF>  ~ 

-  -  •|<DIDUIU>], 

Where  all  expectation  values  are  understood  to  be  taken 
with  respect  to  the  Neel  state* 


2e  The  N-body  Divergence  Problem 


When  we  calculate  explicitly  the  value  of  the 

brackets^  the  brackets  which  contain  a  single  D  operator 

and  a  single  U  operator  are  found  to  be  proportional  to  N, 

the  total  number  of  lattice  sites  in  the  crystal.  Brackets 

.with  more  than  a  single  D  and  U  are  of  a  higher  order  of  N. 

Thus^  all  terms  of  the  which  are  products  of  two  or  more 

2 

brackets  are  at  least  of  the  order  of  N  ^  and  do  not  repre¬ 
sent  a  physical  contribution  to  the  The  argument  given 

for  this  conclusion  is  that  if  E^  then  E  ^  and 

the  energy  per  lattice  slte^  E/Nj  becomes  unbounded  as  the 
number  of  lattice  sites  increases^  l.e.^  llmj^ ^^E/N  — >eo . 
This  is  also  true  for  higher  powers  of  N.  However,  these 
non-physical  terms  are  only  apparent,  for  contained  impli¬ 
citly  in  other  terms  of  the  same  order  we  find  the  negative 
of  these  terms. 

Before  v;e  examine  this  cancellation  for  all  six  or¬ 
ders  of  Eq,  (32),  we  investigate  this  cancellation  in  a 
related  problem,  i.e,,  that  of  a  perturbation  Hamiltonian 
consisting  solely  of  the  operators  D  and  U.  In  other 
words,  we  assume  1=0.  Our  energy  terms  are  then 

E[  =  0, 

E^  =  -  ^DU>, 


(33) 
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E^  =G, 

and 

Eg  =^[-2<PU>^  +  |<D^U^XPU>  -  ^P^UPU^>  -  ^D^u2>]. 

We  proceed  to  calculate  these  terms  explicitly  in  order  to 
show  the  cancellation  of  multiple  bracket  terms. 

In  Eg  we  find  that  <DU>  2  <Q>,  where  Q  =  (P»U),  the 
commutator  of  D  and  U.  For  <Q>  we  find 

<Q>  -  mr.al  (34) 

^0 

In  Ej|^  it  is  found  that  <P^U^>  =  2<Q?>^  +  <P>,  where 
E  =  [D(Q,U)].  a  calculation  of  <F>  yields 

<J>  =  J-a  .a)  X  [6(S^^)2^  -  4(SXq)  +  13. 

°  °  (35) 

In  Eg  we  find  that  <P^U^>  =  6<Q>^  4-  9<QXF>  +  2<R'*'R> 
+  <V>,  where  R  =  {Q#U)  and  V  *■  [D(P,U)],  A  calculation  of 


<R  R>  and  <V>  yields 


<R+R>  = 


V. 


0 


(36) 


0 


?(^^n  +3^.-3) 

+  4(S!rQ)2 - 0 - gO - SCS^Tq)  +  1], 


and 


Al.-a)^!>0',6„,'2  ,  ,4(15»'o"-'‘5J'o+*0) 


<V>  =  X  [4(s^o) 

Or 


'0 


-  120(sro)3^  +  8(s;ro)4-yo^-^36-3) 


0 


(37) 


-38(3^^)  +  7] 


Using  Eqs.  (34) -(37)  In  conjunction  with  Eq.  (33),  the 


energy  correction  through  six  orders  is  found  to  be 


36 


(39) 
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and 


1  4 

-^DI^U>. 

e? 


Since 


IUlO>  =  -  2JUl.0>,  (40) 

these  energy  corrections  are  manifestly  proportional  to  N 
in  view  of  Eq.  (34). 

/  \  4 

For  the  contributions  proportional  to  (1-a)  ,  three 
of  them  have  been  calculated  explicitly,  l.e.. 


4 


<DU>‘ 


}• 


(41) 


<D^UIU> 


12  2 


}> 


and 


-^j4<DUXDI^U>  +  2<DIU>^  -  •|<DIDIU^> 


<D^UI^U> 


-  •|<DIDUIU> 


Using  Eq,  (4o)  and  the  similar  relationships. 


IU^l0>  =  -  2J[2U^I.0>  +  2U^U^,J.0>],  (42) 


where 


U3^Ui,l.o>  =  (|)j2(i-a)2  I 

Aoti 


and 


lVl.O>  =  (-2J)^|4U^I.0>  +  10U3_U^,l0>  + 

(44) 

where 


AoTAt^  A  <r^ 

-  ( V>^) 


we  find  for  ,  and 


=  -  -V^>^ 

A-^ 


(46) 


= 

5 


39 


and 


E, 


a5 


0 


Prom  Eq.  (35) ^  these  energy  corrections  are  also  manifest¬ 
ly  proportional  to  N, 

3o  The  Recovery  of  the  Islng  Model  as 
Zero  Order  Hamiltonian 

Examination  of  Eq,  (39)  Indicates  that  a  very 
regular  series  Is  developing  for  the  terms  proportional  to 

p 

(l-a)  ,  If  we  project  this  regularity  to  orders  higher 
than  the  sixth,  we  obtain  for  the  entire  contribution  to 

p 

the  energy  which  is  proportional  to  (l-a)  , 


Oo 


E 


(l-a)  «  _  i  X  (-  i)''<Dl’^U>  =  - 


'^=0 


(47) 


In  view  of  Eq.  (4o).  Ey  predicting  and  Including  all  terms 

p 

proportional  to  (l-a)  In  the  energy  correction,  we  are  ef¬ 
fectively  recovering  the  Ising  model  as  zero  order  Hamil¬ 
tonian  to  second  order  in  the  perturbation. 

In  order  to  recover  the  .Ising  model  through  fourth 
order,  we  must  be  able  to  determine  the  regi,xlarltles  in  the 
energy  corrections  which  are  proportional  to  (l-a)  ,  which 


we  may  represent  by  the  expression 


nU  ^ 


Assuming  that  terms  proportional  to  powers  of  N  greater 
than  one  are  always  cancelled  out  of  our  energy  correc¬ 
tion  In  each  order  n,  as  they  were  for  n  =  4,  and  6, 
the  physically  Important  part  of  the  energy  corrections  Is 

entirely  contained  In  the  single  bracketed  terms.  These 

/  x4 

terras  have  the  general  form  If  proportional  to  (1-a) 


p,  m,  q  Integer.  (48) 


This  Is  shown  explicitly  below.  It  Is  only  necessary  to 
reduce  the  terms  of  Eq.  (48),  l.e.,  remove  contributions 
proportional  to  powers  of  N  greater  than  one,  to  obtain 
the  non-vanlshlng  contribution  to  the  perturbation  correc¬ 
tion.  If  we  Indicate  a  reduced  term  by  a  subscript  R,  then 
the  entire  contribution  to  the  energy  series  which  Is  pro¬ 
portional  to  (l-a)^  may  be  written 


E 


(1-a)' 


2A' 


>o  ,  p+q-Hn  ,  r^  « 


p,q,m=0 


(^9) 


If  we  let  m  +  p  +  q  =  0,  we  get  the  single  bi>acketed  term 
In  Letting  p  +  ra  +  q  =  1  yields  three  single 
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bracketed  terms^  two  of  which  are  equivalent  because  of  the 
Hermltian  property  of  the  brackets,  l.e,, 

<D^UIU>  =  <D^UIU>‘^  =  «  <DIDU^>.  (50) 

Thus  p  +  q  +  m  =  1  yields  the  single  bracketed  terms  In 
(1-a)^ 

’  ,  Putting  p  +  q  +  m  =  2  yields  the  single  bracketed 
terms  In  provided  we  use  the  Hermltian  properties 

of  the  Individual  brackets. 

With  Eq,  (4o),  together  with  the  relationship 

l”^U^10>  =  (-2J)®-[2™U^t0>  +  2(3“-2’^)U^U3_,  t.0> 

+  (4ni_2.3m^  2"^)U3,2,U^,2l.O>j,  (51) 

we  find  for 


A  similar  treatment  of  the  terms  In  Eq.  (32)  which 
are  proportlona.1  to  (l-a)  may  be  undertaken.  The  expres¬ 
sion  equivalent  to  Eq.  (49)  for  these  terms,  subject  to 
the  assumptions  made  in  obtaining  Eq.  (49),  Is 
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sd-a)®  .  -  1  ‘f  ( 

nijp^  s=0 

X  [<diPdi^di“ui®ui‘^u>p^ 

+  3”^'‘‘^<DI^DI^Ul'"DI®Ul'^U>p].  (53) 

In  Appendix  A  we  give  the  reduced  form  for  Eq.  (53) ^ 
and  perform  the  sums.  The  resulting  expression  for  the 
general  case  is  so  unwieldy  that  It  was  determined  that 
the  additional  correction  obtained  from  calculating  Eq. 
(53)  was  not  worth  the  effort.  It  did,  however^  seem  im¬ 
portant  to  carry  out  the  calculation  of  Eq.  (53)  for  the 
spin  one-half  coupled  linear  chain,  which  is  done  In  Appen 
dlx  A.  The  numerical  result  Is  the  same  as  that  of  a  per- 
turbatlon  series  developed  by  Walker  ^  for  this  particular 
case,  thus  confinning  our  development  of  a  perturbation 
series  with  the  Ising  model  as  zero  order  Hamiltonian. 

C,  The  Perturbation  Series  for  the  Energy 

1.  The  Second  and  Fourth  Order  Corrections 
for  the  General  Cases 

Eor  arbitrary  spin  and  for  the  linear  chain,  plane 
quadratic,  and  simple  cubic  lattices,  the  series  expansion 

O 

In  powers  of  (1-a)  for  the  energy  Is 


^  m+p+q+r+s  ^ 

gr+s  ^ 


E  =  Eq[1  +  C2(l-a)^  +  Ci^(l-a)^  +  Cg(l-a)^  +.  ...](54) 


where 


C 


2 


-  I" 


C4  = 


l6(sr^)[2isrQ)  -  ipmsY^)  -  dUisY^)  - 1] 


K  -1 

W  J. 


X  [192(S^q)5I^  -  l6(S<JrQ) 


46  +  15Yq  -  15 


+  3^0-3  ?22^^o^  +  9i^o-9 

+  32(S>q)3--0 - ^-Q_  -  4(sro)^—  °  ° 


0 


+  24(3^^)  -  3]. 


In  Table  II  are  listed  the  values  of  and  C4  for  the 
three  lattices^  where  the  spin  Is  varied  from  S  =  1/2  to 
3=2,  We  also  give  the  noinnallzed  (to  the  Islng  ground 
state  energy)  energy  for  the  isotropic  case,  a  =  0.  These 
are  necessarily  positive,  and  thus  the  larger  the  numeri¬ 
cal  result,  the  lower  the  corresponding  ground  state  ener¬ 
gy.  For  comparison  we  give  the  results  obtained  by  the 
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linked  spin  cluster  method  of  Davis  for  the  equivalent 
of  Cg  and  C4,  and  for  the  normalized  energy  through  terms 
of  the  order  of  (l-a)^  (Linked  Cluster  I  in  the  table). 

We  also  repeat  the  final  normalized  energy  reported  for 
the  linked  spin  cluster  method,  which  includes  part  of  the 


contribution  from  terms  of  the  order  of  (l-a)  (Linked 
Cluster  II  in  the  Table).  Thus  the  difference  between  the 
final  reported  nomalized  energy  and  the  normalized  energy 
through  terms  of  the  order  of  (l-a)  for  the  linked  spin 
cluster  method  (Linked  Cluster  II  and  I,  respectively)  lies 
with  the  partial  calculation  of  Cg  performed  by  him. 

For  the  linear  chain  with  S  =  1/2,  we  obtain  the  quite 
respectable  value  of  1.7500  for  the  normalized  energy.  The 
exact  energy  eigenvalue  to  five  figures  Is  1.7726."^  This 
may  be  compared  with  the  spin  wave  ground  state  energy  of 
1.7268^^  and  the  linked  spin  cluster  result  of  1.7363. 

Upon  Investigating  Table  II,  we  find  this  work  to  yield 
lower  eigenvalues  than  that  of  the  linked  spin  cluster 
method  for  the  cases  in  which  the  spin,  S,  and/or  the  di¬ 
mensionality  of  the  lattice  Is  small.  The  partial  calcu¬ 
lation  of  C^  In  Linked  Cluster  II  appears  to  be  sufficient 
for  the  two  dimensional  lattice  when  S  >  3/2,  and  for  the 
three  dimensional  lattice  when  S  >  1.  For  the  one  dimen¬ 
sional  lattice,  the  partial  calculation  of  C^j^  does  not 
seem  to  be  adequate. 

2.  The  Energy  Series  In  the  Spin  Wave  Theory 
as  a  Function  of  Islng  Anisotropy 

Turning  to  an  examination  of  the  spin  wave  theory 

of  the  antiferromagnetic  ground  state^^  and  Its  first  or- 

^2 

der  correction  given  by  Oguchl,  we  determine  Its  rela¬ 
tionship  to  the  present  work.  The  spin  wave  theory  has 
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proved  to  provide  the  most  successful  description  of  a 
Heisenberg  antlferromagnet  at  temperatures  T  »  Tjj  ~ 
Neel  temperature)  short  of  an  exact  treatment  of  the  model, 
and  consequently  one  expects  it  to  give  an  adequate  de¬ 
scription  of  the  ground  state.  The  spin  wave  Hamiltonian 
Is  obtained  from  Eq.  (6)  by  replacing  the  angular  momen¬ 
tum  operators  with  equivalent  operators,  written  In  terms 
of  the  operators  applicable  to  the  problem  of  the  linear 
harmonic  oscillator,  i.e.. 


Sj  =  S  -  n^,  (55) 

Sj  =  (2S)^/^(1  -  nj/2S)^/“aj, 

Sj  «  (2S)^/^aJ(l  -  nj/2S)^/^, 

Sj  =  -  S  +  nj^, 

S+  =  (2S)^/2bJ(l  -  nj^/2S)V2, 

Sj^  -  (2S)^/2(1  -  nj^/2S)^/Sj, 


"  W’ 


where 


Substitution  of  the  operators  of  Eq,  (55)  into  Eq.  (6) 
yields 


(56) 


+  2JS(l-a)  I  (1  -  ny2S)^/^a^(l  -  nj^/2S)^/^b 


<Jk> 


"k 


+  2JS(l-a)  I  a*(l  -  n^/2S)^/^bJ(l  -  n^/2S)^^^ 


<Jk> 


V 


2J  Z 
<jk>  ^ 


which,  to  order  1/S  In  the  binomial  expansion  of  the  op¬ 
erators  (1  -  nj/2S)^^^,  etc.,  equals 


* 


+  2JS(l-a)  Z  (a.b.  4a>r) 
<Jk>  ^  ^  ^  ^ 


(57) 
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-  2J  I  n.n,  . 

<Jk>  ^ 

The  first  three  terms  of  Eq.  (57)  constitute  the  spin 
wave  approximation,  and  the  last  two  teirnis  constitute  the 
perturbation  Hamiltonian  considered  by  Oguchi.  If  the 
first  terra  of  Eq.  (57)  is  considered  to  be  of  order  unity, 
then  the  spin  wave  approximation  is  correct  to  order  1/S, 
and  all  of  Eq.  (57)  is  correct  to  order  1/S  . 

We  repeat  the  spin  wave  calculation,  and  the  Oguchl 
correction,  to  the  ground  state  because  of  the  insertion 
of  the  anisotropy  parameter  (l-a)  in  Eq.  (57).  It  will 
subsequently  be  shown  that  when  a  =  0,  we  o]3taln  the  re¬ 
sults  reported  by  Anderson  and  Oguchi.  For  the  purpose  of 
extracting  the  long  range  order  parameter  from  our  expres¬ 
sion  in  later  sections,  we  add  to  Eq.  (57)  the  additional 
anisotropy 

^  "k'- 

(58) 

Eq.  (58)  represents  the  anisotropy  provided  by  an  effec¬ 
tive  field  which  is  directed  oppositely  on  each  sublattice, 
but  it  is  Inserted  here  simply  because  it  allows  the  sub¬ 
lattice  magnetization  to  be  determined  from  the  energy 

PQ 

series  by  means  of  the  Feynman  theorem.  ^ 

The  diagonalization  of  the  spin  wave  Hamiltonian  is 
accomplished  by  means  of  two  successive  canonical  trans- 


formations.  The  first  transformation  is  to  the  operators 
belonging  to  the  reciprocal  lattice. 


(M-J)aj,  (59) 

exp  (-IA-J)a* 

~  (jj)  ^  ^  ^ 

V=  e^P  (5/^-k)b*. 

With  these  transformations;,  Eq.  (57)  plus  Eq,  (58)  be¬ 
comes 

n/ 

H  =  -  2JS^(~^)  +  2g^SH^(|)  (60) 

•jf’ 

+  2JS^q(1  +  2JSK  )  t  L  3 

^  A 

+  2JS(l-a)  I  )^(a^b^+a*b*) 

-  2 J ( 1  -a ) (|) ^ [ 5 ( - /^g- A3+^4 ) 


+  6(-Aj^+A2-A3+Aj,)?^^a*^b*^a*^a^^ 

+  S  {-Ai+A2+-l3-Aj,)A^^^,^,^.^^a*^l,^b*^b^^] 


Using  the  transformation  to  spin  wave  creation  and  destruc¬ 
tion  opera  tors  j, 


cosh©  slnh0j^5  (6l) 

=  -  (A^  slnh©  +  A  cosh  ©^, 
a^  =  0(^  cosh0p^^  -  slnh0;^, 

^  ‘X* 

b^  =  -  slnh0p(^  +  cosh©^ 

with 

If^il-a)  g^. 

^anh^®A  “  ^  D  i  -  1  +  2JS^  " 

we  can  diagonalize  the  spin  wave  Hamiltonian  [the  first 
four  terms  of  Eq,  (6o)],  resulting  In 

=  2Js2(^)  +  2g^^(f)S 


(62) 
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A  ®o 


2JS^oIiD  "  C®‘ 
A 


2  2rx^(i”a)2  1/2 


Y  2 
^0 


-] 


ft 


^  w 

For  the  ground  state^,  <^^0\>  =  0  and  =  0,  and 

A 


«■ 


if  also  H«  =  0  we  get  for  the  energy  of  the  ground  state 


N3/ 


2JS^(-2^)  -  2JSrQ^|l 

A  ^ 


[1 


^^‘^(l-a)'^  1/2) 


r 

(63) 


which  is  Anderson's  result  when  a  =  0, 

If  we  compute  the  first  order  perturbation  correction 
to  Eq,  (63)  using  the  remaining  terms  of  Eq„  (60)  as  a 
perturbation  Hamiltonian^,  we  obtain  the  equivalent  of  the 
ground  state  energy  correction  obtained  by  Oguchl^  but 
with  the  insertion  of  the  anisotropy  parameter  (l-a). 
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which  yields  Oguchl's  correction  when  a  =  0. 

Table  III  gives  the  eigenvalues  E  +  AE  for  the  cases 
Investigated  In  Table  II.  In  most  cases,  the  spin  wave 
theory  plus  the  Oguchl  correction  yields  eigenvalues  which 
are  lower  than  those  yielded  by  our  method.  However,  the 
correction  for  the  linear  chain  with  S  =  1/2  badly  over¬ 
shoots  the  exact  eigenvalue,  while  for  the  same  lattice 
with  S  =  1  our  method  gives  an  eigenvalue  which  lies  lower 
than  that  of  the  corrected  spin  wave  theory.  This  oscilla¬ 
tory  behavior  of  the  corrected  spin  wave  theory  seems  to 
be  a  property  of  the  smallness  of  the  spin,  for  a  similar 
behavior  may  be  seen  In  the  plane  quadradlc  and  simple  cu¬ 
bic  lattices.  Of  course  there  are  no  exact  values  with 
which  to  compare  in  the  two  and  three  dimensional  cases, 
and  the  oscillatory  behavior  is  Inferred  from  the  behavior 


of  the  corrected  spin  wave  eigenvalues  relative  to  those 
of  our  method. 

3.  Comparison  of  the  Energy  Series  from  the 

Spin  Wave  Theory  with  that  of  the  Perturbation 
Theory 

The  relationship  between  the  spin  wave  approach 
and  the  modified  perturbation  approach  may  be  determined 
in  a  more  elaborate  manner  than  by  a  simple  comparison 
of  eigenvalues.  This  is  possible  because  of  the  presence 

p 

of  the  two  expansion  parameters^  1/S^q  and  (1-a)  ^  common 
to  the  energy  series  of  both  methods.  Ihe  entire  contribu 
tlon  to  the  normalized  energy  eigenvalue  of  the  antiferro¬ 
magnetic  ground  state  resulting  from  a  continuation  of 
either  of  the  above  methods  may  be  expressed  in  the  follow 
Ing  general  form. 


OO 

E/E  »  1  +  I  a  (l-a)^(l/S^Q)“,  (65) 

n,  m=l 

Examination  of  Eq.  (54)  and  Eq.  (65)  reveal  the  fol¬ 
lowing  relationships. 
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and  examination  of  Eq„  (63) ^  Eq.  (64)^  and  Eq,  (65)  yields 


Spin  Wave  Correction  = 

iiSQ 


(67a) 


Oguchl  Correction  or  j 
First  Order  Perturbation L 
Correction  to  the  Spin  [ 
Wave  Theory 


p  ^  a„  p  ( 1  "a ) 

{sSgf  n-l 


2n 


(67b) 


Prom  Eq.,  (54),  we  can  develop  the  first  few  terms  of 
Eq.  (66)  explicitly,  and  from  Eq,  (63)  and  Eq,  (64)  we  can 
develop  the  first  few  terms  of  Eqs,  (67a,b)  explicitly. 

The  terms  which  are  common  to  Eq.  (66)  and  Eqs,  (67a, b) 
are 


an  =  1/2,  (68) 
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Eq.  (65)  may  be  expressed  In  a  tabular  form.  In  which 


o 

the  rows  are  labeled  according  to  powers  of  (l-a)  ^  and  the 
columns  by  powers  of  I/S'Sq. 
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Spin  wave  theory  and  Its  successive  corrections  cor- 
respond  to  adding  the  terms  In  the  next  column  to  the  pre¬ 
vious  results.  Beginning  with  the  Islng  model  for  a  zero 
order  Hamiltonian  and  computing  succeeding  orders  of  the 
perturbation  corresponds  to  taking  Into  account  the  terms 
of  successive  rows, 

4,  A  Hybrid  Theory 

It  now  becomes  possible  to  take  advantage  of  all 
of  the  terms  In  Eq.  (69)  In  both  the  first  two  columns  and 
the  first  two  rows  to  compute  the  ground  state  energy  ei¬ 
genvalue.  This  Is  done  simply  by  adding  together  C2j^^ 

spin  wave  correction^,  and  first  order  (Oguchl)  correction. 
The  twice  counted  terms^,  aj^2f  ^21^  ^22^  then  sub¬ 

tracted  from  the  resulting  sum.  These  results  are 


presented  in  Table  III  under  the  heading  **%brid  Method, " 
The  hybrid  method  yields  an  eigenvalue  lying  lower  than 
either  the  corrected  spin  wave  or  the  method  of  this  dls- 
sertatlon^,  except  for  the  spin  one-half  cases  and  the  sim¬ 
ple  cubic  spin  two  case. 

The  corrected  spin  wave  theory  may  be  expected  to  ap¬ 
proach  the  exact  eigenvalue  for  large  spin  and  number  of 
nearest  neighbors.  For  the  spin  one-half  cases,  however, 
the  apparently  "better"  (l.e.,  lower)  eigenvalues  produced 
by  the  corrected  spin  wave  theory  may  be  misleading.  In 
the  spin  one “half  linear  chain,  the  lower  eigenvalue  given 
by  the  corrected  spin  wave  theory  must  be  raised  by  addi¬ 
tional  corrections,  because  it  lies  a  great  deal  lower  than 
the  exact  eigenvalue.  In  this  particular  case,  the  hybrid 
theory  yields  an  eigenvalue  which  is  less  than  in  error. 
It  seems  likely  that  because  of  the  behavior  of  the  cor¬ 
rected  spin  wave  eigenvalues  which  was  pointed  out  at  the 
end  of  Part  (2),  Section  C  of  this  chapter,  the  hybrid 
method  gives  the  best  eigenvalues,  except  as  the  spin  gets 
relatively  large. 

We  now  turn  to  the  consideration  of  the  short  and 
long  range  order  parameters  for  the  antiferromagnetic 
ground  state. 


D.  Short  and  Long  Range  Order 
1.  The  Feynman  Theorem 

There  is  a  theorem  due  to  Peynman^^  which  we  now 
need  in  order  to  detemiine  a  perturbation  series  for  the 
short  and  long  range  order  parameters  from  the  series  for 
the  energy.  Brief ly^  the  theorem  is  thlss  Given  the 
eigenvalue  equation 


H(A)Cj;(A)  -  U(A)(p(A).  (70) 

where  H(/\)  is  an  operator  expressed  as  a  function  of  a 
parameter  At  then 


SUfA) 

3A 


(71) 


For  a  proof  of  Eq,  (71 )>  the  reader  is  referred  to  the 
article  by  Feynman. 

2.  The  Short  Range  Order  in  the  Modified 
Perturbation  Theory 

Once  the  short  range  order  parameter  has  been 
defined,  the  Feynman  theorem  may  be  used  to  generate  the 
corresponding  series  from  the  energy  series.  However,  it 
is  necessary  that  the  Hamiltonian  be  expressed  as  an  ex¬ 
plicit  function  of  the  short  range  order  operator.  The 
natural  definition  of  the  short  range  order  operator  is 
one  which  measures  local  spin  correlation  along  some  pre- 


f erred  axlSj,  i,e,j 


1  T  qZqZ 

^  i-S) 


>]op  is  the  average  spin  correlation  along  the  axis  of  an¬ 
isotropy  between  spin  systems  on  neighboring  lattic  sites. 
It  is  noimiallzed  so  that  for  perfect  ferromagnetic  ordering, 
<)^p>  »  1^  for  a  system  of  conpletely  independent  spin  sys¬ 
tems,  <^Qp>  “  0,  and  for  perfect  antiferromagnetic  order¬ 
ing,  <7|qp>  =  -1. 

In  terms  of  the  Hamiltonian  of  Eq,  (6),  and  thus  also 
of  Eq.  (21),  the  short  range  order  operator  takes  the  form 


7inr[i 

_ y  \ 


Application  of  the  Feynman  theorem  to  Eq.  (5^)  then  yields 


bp  true  ground  state 


213*^  (-^) 


5fjr-[l  +  (l-a)^][Eo[l  +  Og(l-a)' 


+  CkCi's)  +  , . . ] 


-  1  +  Cgd-a)^  +  SG^Cl-a)^  +  5C5(l-a)^  + 


Values  of  ^  in  the  Isotropic  liml.t  (aO)j,  determined,  from 
Eq„  {Y^)s  Sire  given  in  Table  IV,  It  is  interesting  to  com¬ 
pare  the  approximate  value  of  determined  from  Eq,  (7^)s 

with  the  exact  value  of  in  the  Isotropic  limits,  deter- 
l4 

mined  by  Orbach^  for  the  spin  one-half  coupled  linear 
chain.  These  are^  respectively ![=  -0,7500  and 
-0.5909.  ThuSj  Eq,  (74)  does  not  represent  the  short  range 
order  with  nearly  the  accuracy  with  which  Eq.  (54)  repre¬ 
sents  the  energy  to  this  order  in  the  perturbation.  In 
Appendix  A  we  determine  Cg  to  be  zero  for  the  spin  one- 
half  coupled  linear  chain^  so  that  the  value  of  is  not 
further  modified  until  the  eighth  order  of  the  perturbation. 
In  Section  E  of  this  chapter,  we  will  discuss  further  the 
properties  of  the  short  range  order  series  for  this  par¬ 
ticular  case  in  conjunction  with  a  discussion  on  the  con¬ 
vergence  properties  of  the  energy,  short  range  order,  and 
long  range  order  series  for  the  general  case, 

3.  The  Short  Range  Order  in  the  Corrected 
Spin  Wave  Theoiry 

Application  of  Eq,  (73)  to  the  spin  wave  perturba¬ 
tion  expansion  of  the  energy,  Eq.  (63)  plus  Eq„  (64), 
yields  the  short  range  order  as  determined  by  the  corrected 
spin  wave  theoi*y. 


1  S2r^ 
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7^  =  -  1  +  G. 
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Values  of  in  the  isotropic  limit  (a  =  O)^  detei^nlned  from 
E<1«  (75)  through  Gg,  are  given  in  Table  IV.  However^,  no 
values  are  given  for  in  this  approximation  for  the  linear 
chain.  This  is  because  Eq.  (75)  is  indeterminate  for  the 
one  dimensional  lattice.  Both  G^^  and  Gg  are  not  finite  for 


this  casGj  and  thus  this  indeterminacy  is  independent  of 
the  spin.  This  behavior  is  definitely  a  property  of  the 
approximation^  and  not  of  the  short  range  order  parameter, 
since  from  Orbach's  work  we  know  the  short  range  order 
in  the  Isotropic  limit  of  the  spin  one “half  coupled  linear 
chain  to  ,be  finite.  We  will  find  this  behavior  to  be  true 
also  of  the  long  range  order  for  this  lattice,  when  it  is 
determined  by  means  of  the  spin  wave  approximation.  It  is 
most  likely  due  to  the  convergence  properties  of  the  series 
for  these  cases,  i.e.,  when  the  terms  in  the  series  for  the 
short  and  long  range  order  are  expressed  in  a  form  similar 
to  that  of  Eq.  (69)  for  the  energy,  convergence  depends  up¬ 
on  the  order  in  which  -che  terms  are  added  up.  If  they  are 
added  up  by  rows  (modified  perturbation  theory),  conver¬ 
gence  is  evident,  and  if  they  are  added  up  by  colvumis 
(spin  wave  approximation),  the  series  diverges. 

In  Table  IV,  we  also  Include  the  short  range  order 
parameter  as  determined  by  the  hybrid  method  of  Section  5j 
Part  (4)  of  this  chapter. 

4.  Comparison  of  the  Methods  for  the 

Determination  of  the  Short  Range  Order 

In  view  of  the  fact  that  the  corrected  spin  wave 
theory  is  Inadequate  for  the  determination  of  the  short 
range  order  of  the  only  case  for  which  the  exact  value  is 
known,  no  absolute  comparison  of  the  methods  is  possible 
at  present.  Since  the  spin  wave  approximation  represents 
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an  expansion  in  powers  of  I/S^q^  one  expects  it  to  be  the 
best  method  when  the  spin  and  dimensionality  of  the  lattice 
are  large.  However as  has  already  been  pointed  out^  the 
short  range  order  series  diverges  for  the  linear  chain  no 
matter  how  small  the  expansion  parameter^  and  thus  one  can¬ 
not  fully  depend  upon  this  criterion.  We  conclude  that, 
for  the  short  range  order  parameter,  the  modified  perturba¬ 
tion  theory  gives  the  more  meaningful  prediction.  This 
conclusion  will  be  borne  out  in  the  discussion  of  the  rel¬ 
ative  convergence  of  the  series  generated  by  the  two  meth¬ 
ods  which  appears  in  Section  E  of  this  chapter, 

5.  The  Long  Range  Order  in  the  Modified 
Perturbation  Theory 

The  choice  of  an  operator  to  represent  the  long 
range  order  is  restricted  here  by  the  conditions  necessary 
for  the  application  of  the  Feynman  theorem.  The  natural 
extension  of  the  operator  which  was  used  to  represent  the 
short  range  order  would  be  one  which  measured  the  spin 
correlation  along  the  axis  of  anisotropy,  l„e,, 

5”^oo 

or  one  which  measures  the  general  spin  correlation, 

/^op  “  S(S+iy^j'^j+r 


(77) 
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The  correlation  is  measured  between  spin  systems  which  lie 
at  large  distances  from  each  other^  hence  the  concept  of 
long  range  order.  When  a  wave  function  is  available^  such 
definitions  of  the  long  range  order  operator  as  Eqs„  (76- 
77)  are  possible^^^  but  when  use  is  made  of  the  Feynman 
theorem^  the  condition  that  the  Hamiltonian  must  be  an 
explicit  function  of  the  long  range  order  operator  pre¬ 
vents  the  adoption  of  either  Eq.  (76)  or  Eq,  (77).  Con¬ 
sequently,  we  are  constrained  to  use  the  concept  of  the 
sublattice  magnetization,  which  bears  some  relationship  to 
Eq,  (76). 

The  sublattice  magnetization  operator  may  be  expressed 
as  follows; 


Thus,  for  perfect  order,  <p>  =  ±  1,  depending  upon  the  di¬ 
rection  in  which  the  sublattice  is  magnetized.  For  com¬ 
plete  disorder,  <p>  =  0.  The  concepts  of  order  and  disor¬ 
der  used  in  conjunction  with  Eq,  (78)  are  not  identical 
with  those  used  with  Eq.  (76).  The  correlation  between 
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spins  at  large  distances  can  be  such  as  to  yield  zero  for 
the  expectation  value  of  Eq.  (78);,  while  yielding  a  finite 
value  for  the  expectation  value  of  Eq.  (76) .  In  other 
wordSj  the  spins  at  individual  lattice  sites  may  be  chang¬ 
ing  with  time  in  such  a  way  as  to  wash  out  the  average  2 
component  of  spin,  and  yet  be  in  step  with  the  time  varia¬ 
tion  of  each  other  so  that  the  average  product  is  non¬ 
zero.  If,  however,  the  sublattice  magnetization  is  finite, 
then  the  average  z  component  of  spin  at  individual  lattice 
sites  is  not  zero,  and  the  spins  are  "pinned  down"  to  point 
mainly  in  the  z  direction.  In  this  case,  Eq.  (78)  bears 
the  following  approximate  relationship  to  Eq.  (76), 


We  will  use  this  relationship  later,  when  we  compare  values 
of  P  with  7|^, 
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In  order  to  apply  the  Feynman  theorem^  It  is  necessary 
to  express  the  effective  Hamiltonian^  Eq.  (21)^  as  an  ex¬ 
plicit  function  of  Eq.  (78).  We  repeat  Eq,  (21)^  with  Hq 
replaced  by  Its  eigenvalues  with  respect  to  the  set  of 
states  described  in  Eq.  (18)^ 
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The  difference  In  the  forms  of  Dj^  for  S^q  =  1  and  S^Tq  7^  1 
arises  from  the  fact  that  in  Eq.>  (52)  <^^%2'^1'2^R  ”  ® 
for  =  1,  and  finite  for  SSq  7^  1. 

6.  The  Long  Range  Order  in  the  Corrected 
Spin  VJave  Theory 

Turning  again  to  the  corrected  spin  wave  theory, 
we  determine  the  expectation  value  for  the  long  range  or¬ 
der  operator,  Eq.  (78).  In  view  of  Eq,  (58)  and  the 
Hamiltonian,  Eq.  (60),  application  of  the  Feynman  theorem 
yields 


r  gps 
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where  E  Is  the  sum  of  Eqs,  (63)  und  (64),  The  series  for 


pis  then 
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Where 


Kn  o  \  o  “V2 

%  =  t(l  -  -  1]. 
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For  a  -  Oj  K-^  Is  the  spin  wave  correction  to  the 
normalized  sublattice  magnetization  (l,e.j  our  long  range 
order  parameter),  and  Kg  =  0,  consistent  with  the  result 
of  Oguchl.^^ 

7o  Comparison  of  the  Methods  for  the  Determination 
of  the  Long  Range  Order 

In  Table  V  are  given  the  values  of  Dg  and  for 
the  cases  Investigated  previously.  In  Table  VI  are  given 
the  values  of  ^  for  the  isotropic  case  (a=0)  for  the  same 
cases,  listed  under  the  heading  "this  work,"  Also  given 
in  Table  VI  are  the  values  of  P  predicted  by  the  corrected 
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spin  wave  theory  and  the  hybrid  method  of  Part  (4)^  Sec¬ 
tion  G  of  this  chapter.  No  values  of  p  appear  for  the 
latter  two  methods  as  applied  to  the  linear  chaln^  for 
reasons  similar  to  those  given  for  the  omission  In  Table 
IV  with  respect  to  the  determination  of  7^[cf,  Part  (3)5 
Section  D  of  this  chapter].  This  omission  will  be  dis¬ 
cussed  further  In  Section  E. 

The  total  sublattice  magnetization^,  determined  by 
either  the  modified  perturbation  theory  or  the  corrected 
spin  wave  theory,  may  be  expressed  In  a  form  similar  to 
that  In  which  the  energy  was  expressed  in  Eq,  (69)^.  l.e.. 
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Which  are  determined  from  Eq.  (83)  and  either  Eq„  (84)  or 
Eq»  (85)i  or  alternatively  from  Eq.  (87)0 

In  order  to  determine  the  power  of  the  spin  wave  ap¬ 
proach  versus  the  modified  perturbation  theory^  we  also 
Include  in  Table  VI^  under  the  heading  of  "overlap, "  the 
contribution  to  p  which  arises  solely  from  the  four  terms 
common  to  both  methods,  i.e.. 


overlap  -  1  -  («iii+d21^s|^  "  ^‘^12‘‘‘‘^22^  ^2^  ^2  “ 

°  (90) 

It  is  apparent  from  Table  VI  that  the  overlap  between 
the  two  methods  constivutes  far  and  away  the  largest  con¬ 
tribution  to  the  divergence  of  the  sublattice  magnetiza¬ 
tion  from  unity  in  the  three  dimensional  cases.  Not  quite 
so  strong  a  statement  may  be  made  for  the  two  dimensional 
cases,  and  the  modified  perturbation  method  seems  to  be 
the  only  satisfactory  approach  for  the  one  dimensional 
lattice.  Thus,  both  methods  are  roughly  comparable  in 
three  dimensions,  and  one  could  not  hope  to  obtain  very 
different  results  for  the  ground  state  sublattice  magneti¬ 
zation  by  choosing  one  method  in  preference  to  the  other. 

The  long  range  order  parameter  is  larger  in  absolute 
magnitude  than  the  short  range  order  parameter  in  almost 


all  of  the  cases  of  the  Tables  IV  and  ¥!»  At  first  glance. 
It  might  be  thought  that  ^  should  at  all  times  be  larger 
than  p,  because  of  the  comparative  distances  over  which 
correlation  is  being  measured.  The  fact  that  this  is  not 
reflected  in  the  computed  values  m^y  be  attributed  to  the 
way  in  which  our  long  range  order  parameter  has  been  de¬ 
fined.  Our  definition,  Eq.  (78),  which  is  the  only  one 
available  to  our  theory,  is  not  the  proper  analog  of  our 
short  range  order  parameter,  Eq.  (72),  to  lend  Itself  di¬ 
rectly  to  a  comparison  with  This  was  pointed  out  in 
Part  (4)  of  this  section  when  we  derived  an  approximate 
relation  between  the  correct  analog,  Eq.  (76),  and  Eq. 

(78).  For  comparison  with  the  short  range  order  parameter, 

p 

we  have  included  a  table  of  p  ,  Table  VII,  which  is  based 
on  the  approximate  relation,  Eq.  (79). 

8.  Comparison  of  the  Theoretical  Long  Range 
Order  with  Experiment 

There  are  two  experimental  determinations  of  the 
sublattice  magnetization  in  real  antlferromagnets  in  the 
literature.  They  correspond  to  the  measurement  of  the 
sublattice  magnetization  in  a  body-centered  cubic  spin 
5/2  system, (the  magnetic  lattice  in  MaPg),  and  of  the 
sublattice  magnetization  in  a  simple  cubic  spin  5/2  sys¬ 
tem,  (the  magnetic  lattice  in  KMnP^)* 

The  body-centered  lattice,  while  of  the  two  sublat¬ 
tice  stJTUcture,  has  not  been  discussed  heretofore  because 


of  Its  essentially  different  nearest  neighbor  structure. 

In  the  three  lattices  which  have  been  considered^  i.e., 
linear  chains,  plane  quadratic,  and  simple  cubic,  the  posi¬ 
tions  of  the  nearest  neighbors  define  an  orthogonal  coor¬ 
dinate  system,  and  hence  the  6  are  linearly  Independent. 
This  is  not  true  of  the  S  corresponding  to  the  body- 
centered  lattice,  and  hence  greater  care  must  be  taken  in 
determining  the  valties  of  the  sums  over  the  reciprocal 
lattice  which  lead  from  Eqs.  (47)  and  (52)  to  the  values 
of  the  Cg  and  of  Eq.  (54).  QSius,  substitution  of 
JTq  =  8  in  Eq.  (54)  does  not  lead  to  the  values  of  Cg  and 
Gj^  which  are  characteristic  of  the  body-centered  lattice. 

As  a  result,  we  omit  discussion  of  this  case  here,  except 
to  say  that  there  is  an  apparent  discrepancy  between  the 
experimental  result,  ((>=1,00±0,02),  and  the  spin  wave  pre¬ 
diction,  ((>=0.976) ,  This  same  type  of  discrepancy  appears 
in  the  simple  cubic  case,  (KMnF^)^  and  here  direct  com¬ 
parison  with  the  theory  as  presented  in  this  dissertation 
can  be  made. 

The  experimentally  measured  sublattice  magnetization 
in  KMnP2  has  been  determined  to  be  p=  O.998  ±  0,015. 

From  Table  YI,  we  obtain  for  the  results  of  this  disser¬ 
tation,  p=  0.977j  and  for  the  spin  wave  theory,  ()=  O.969, 
The  results  of  the  spin  wave  theory  are  clearly  in  dis¬ 
agreement  with  the  experimental  results,  while  the  results 
of  this  work  lie  Just  outside  the  possible  experimental 


error.  With  the  calculation  of  additional  terms  of  the 
perturbation  series,  Eq.  (82),  it  might  be  hoped  ^  would 
converge  to  a  limit  lying  within  the  interval  indicated 
by  experiment.  That  such  is  probably  not  the  case  is 
shown  by  the  following  argument. 

We  consider  the  contribution  made  to  the  terras  Dg  and 
appearing  in  the  perturbation  series  for  the  sublat¬ 
tice  magnetization,  i.e.,  Eq.  (5^)j  by  the  teimis  of  the 
overlap,  Eq,  (90).  Thus  we  have  for  Dg  and  in  this  ap¬ 
proximation,  using  the  results  of  Eq.  (90),  and  putting 
S  =  5/2,  ^0  =  6, 


TV  ^  1  ^ 

^2  ~ 


^12 


=  0.0178 


(91) 


and 


D4 


1  .. 

S^21 


=  0.0G50. 


Comparison  of  these  approximate  values  with  the  exact 
values  of  Table  V  show  Dg  and  to  be  given  correctly  by 
the  overlap  terms  to  within  2^.  If  we  assume  this  rela¬ 
tionship  to  hold  for  arbitrary  D^,  then  an  upper  limit  to 
the  deviation  from  perfect  order  lies  within  2^  of  the 
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corrected  spin  wave  deviation^  and  a  lower  limit  to  the 
deviation  is  given  by  the  present  theory  to  fourth  or- 
der.  Thus, 


0.968  < 


p  <  0.977. 


(92) 


Of  course,  Eq,  (92)  is  by  no  means  rigorous,  but  it  seems 
to  reflect  the  general  trend  of  the  perturbation  series 
for  the  sublattice  magnetization. 

The  disagreement  of  the  sublattice  magnetization 
predictions  of  both  the  corrected  spin  wave  theory  and 
the  perturbation  theory  with  those  of  experiment  thus  is 
not  resolved.  Considerations  of  canting  in  KMnP^,  as 
well  as  of  anisotropy,  fall  to  account  for  the  difference. 
Canting  has  been  shown  to  have  a  negligible  effect  on  the 
spin  wave  spectrum  for  the  degree  of  canting  present  in 
KMnP^.^^  fhe  anisotropy  field,  Hj^,  in  KMnP^  is  of  the 
order  of  4  oe,  while  the  exchange  field,  Hg,  is  of  the 
order  of  1.6  x  10^  oe.^"^  The  relationship  between  the  an¬ 
isotropy  parameter  a,  H^,  and  Hg  for  the  sublattice  mag¬ 
netization  is  given  by 


a 


+  H 


"E 


‘A 


(93) 


Eq.  (93)  is  based  upon  the  way  in  which  a,  and  Hg  en¬ 
ter  into  the  spin  wave  expression  for  the  sublattice 


magnetization.  Thus  a  Is  of  the  order  of  10“°  for  KMnP2^ 
clearly  too  small  to  affect  our  determination  of  p. 

Experiments  involving  the  determination  of  the  sub- 
lattice  magnetization  for  spin  one-half  systems,  e.g,., 
CuCl2*2H20,  where  the  theoretical  deviation  is  much  larg¬ 
er,  may  make  more  clear  how  to  correlate  theory  and  exper- 

I 

iment, 

E,  Convergence  Properties  of  the  Ground  State 
Parameter  Series  Approximations 

To  begin  with,  almost  nothing  can  be  said  about  the 
convergence  of  series  for  which  the  form  of  the  general 
term  is  not  known.  Of  the  antiferromagnetic  ground  state 
parameter  series,  the  form  of  the  general  term  is  known 
only  for  the  energy  and  short  range  order  series  corres¬ 
ponding  to  the  linear  chain  with  spin  one-half, In  this 
case,  both  series  are  very  noisy,  l.e,,  have  no  apparent 
pattern  for  the  sign  of  succeeding  terms,  and  the  conver¬ 
gence  properties  are  not  readily  assessable.  For  the 
perturbation  series  developed  in  this  dissertation,  the 
discussion  of  the  convergence  properties  will  have  to  be 
limited  to  comparison  with  the  spin  wave  series,  along 
with  comments  on  the  relative  magnitudes  of  succeeding 


terms 


1.  The  Spin  Wave  Series 


We  list  below  the  series  approximations  obtained 
from  the  corrected  spin  wave  theory  for  the  ground  state 
parameters.  The  subscripts  on  the  parameter  ssonbols  re¬ 
fer  to  the  dimensionality  of  the  lattice^  l.e.,  linear 
chain,  plane  quadratic,  and  simple  cubic, 

^  =  1  +  G.7268(|^)  +  0.1321(^)^  +  ...,  (94) 

°D*1 

=  1  +  0.316(^)  +  0.025(-^)^  +  ..., 

°D=2 

^  =  1  +  0.194(^)  +  0.009('J5)^  +  ..., 

°D-3 

^D=2  “  -  1  +  0.786(i)  -  0.756(^)^  +  ...,  (95) 

=  ~  1  +  0, 312('2g)  ~ 

Pd=2  "  ^  ■  0-394(^)  -  O.OOO(^)^  +  ...,  (96) 

1  ■  0»156('^)  -  O.OOO(^)  + 


•  «  •  ^ 
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2.  The  Perturbation  Series 

For  comparison  with  Eqs.  (94)j.  (95),  and  (96), 
we  pick  the  spin  one-half  and  spin  two  cases  from  Tables 
II  and  V. 


^0  D=1 
S=l/2 


=  1  +  (1-a)^ 


-0.2500  (1-a) 


4 


+ 


•  «  •  ^ 

(97) 


J- 

®0  D=1 
S*=2 


=  1  +  0.l429(l-a)^  +  0.0456(l-a)^  + 


®0  I>=2 
S»l/2 


1  +  0.0333(l-a)^  -  0.0065(l-a)^  + 


J- 

^0  D=»2 
S»2 


1  +  0.0667(l-a)^  +  0.0078(l-a)^  +  ..., 


D=3 

S«l/2 


1  +  0.2000(l-a)^  -  0.0008(l-a)^  + 


"*0  0=3 
S-2 


=  1  +  G.o435(l-a)^  +  G.0037(l-a)^ 


+ 


•  »  d  > 


^]>2  =  -  1  +  0.3333(l-a)^ 

3=1/2 


0.0195(l-a)^  + 

(98) 


7|  ^  0.o667(l-a)^  +  G.0234(l-a)^  + 

3=2 


^D=3  =  -  1  +  0.2000(l-a)^  -  0.0024(l-a)^  +  ..., 

3=1/2 

=  -  1  +  0.o435(l-a)^  +  o.oill(l-a)^  +  .... 

3=2 


^  “  0.2222(l-a)^  -  0.0356(l-a)^ 

3=1/2 


•  •  •  i 

(99) 


Pp^2  =  1  ■  0.0356{l-a)^  -  0.0136(l-a)^  - 
3=2 


pP^2  =  1"  0,1200(l-a)^  -  0.008l(l-a)^  - 
3=1/2 

=  1  -  0.0227(l-a)^  -  0.006l(l-a)^  -  ... 

3=2 


The  apparent  convergence  of  the  perturbation  series  Is 
relatively  faster  than  that  of  the  corresponding  spin 
wave  series  when  the  spin  Is  small.  This  judgment  Is 
made  solely  on  the  basis  of  the  magnitude  of  the  ratio 


formed  from  the  last  two  terms  of  each  series.  The  nota¬ 
ble  exception  to  this  rule  occurs  In  the  case  of  the  sub¬ 
lattice  magnetization  parameter,  where  the  first  order 
correction  to  the  spin  wave  theory  Is  In  all  cases  zero. 
Higher  order  corrections  are  not  necessarily  zero,  how¬ 
ever,  and  If  they  were  available,  the  relative  convergence 
could  -be  discussed  for  this  parameter  In  terms  of  them. 


IV.  PERTURBATION  THEORY  OP  THE  HEISENBERG 
ANTIPERROMAGNET  AT  LOW  TEMPERATURES 

A.  The  Partition  Punctlon 

We  now  turn  to  the  consideration  of  the  partition 
function  for  a  Heisenberg  antiferromagnet.  Kubo^^  has 
given  a  perturbation  expansion  of  the  partition  function 
which  is  particularly  suited  to  the  extension  of  the  pre¬ 
vious  treatment  of  the  ground  state  to  include  excited 
states  as  well.  According  to  Kubo, 


.  c»«> 

exp  (-/H)  .  _^(h. 


m=0 


1  E  - 


«0 


)  dE, 

(100) 


where  H  = 


«G 


+  E 


1' 


1.  Zero  Order  and  Perturbation  Hamiltonians 

Since  is  an  operator,  we  need  a  definition  for 
the  operator  1/(E-I^).  We  shall  be  concerned  with  the 
trace  of  Eq,  (lOO)  with  respect  to  a  complete  set  of  ei¬ 
genfunctions  of  for  which  the  definition  of  1/(E-I^) 
is  immediate.  If  we  identify  the  zero  order  and  perturba¬ 
tion  Hamiltonians  with  the  division  of  the  Hamiltonian 
used  in  the  ground  state  perturbation  calculation  of  Chap¬ 
ter  III,  then  we  have  the  division  indicated  in  Eq.  (21). 
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We  may  think  of  Hq  and  as  expressed  in  lattice 
space  or  in  reciprocal  lattice  space^  loe«.  In  the  form 
of  Eq.  (2)j  Chapter  II,  or  of  Eq.  (22),  Chapter  III,  re¬ 
spectively.  For  convenience  In  defining  a  complete  set 
of  eigenfunctions  for  Hq,  we  use  the  Hamiltonian  In  the 
form  of  Eq,  (2).  For  calculatlonal  purposes,  however,  we 
will  Insert  the  transformation  to  reciprocal  lattice  space. 
Thus  the  eigenfunctions  of  Hq  all  take  the  form 

^(n,n^n2)  =  |.n,n^n2>  (lOl) 

where  jn  nj^n2>  represents  a  normalized  eigenstate  of  Hq 
In  which  there  are  n  deviations  from  the  Neel  state,  n^^ 
of  them  on  sublattice  #1,  and  ng  of  them  on  sublattice  #2. 
Thus  n^  +  Hg  =  n.  We  do  not  here  explicitly  differentiate 
among  the  many  different  configurations  of  the  n^^  and  ng 
deviations,  although  this  difference  will  figure  Impor¬ 
tantly  In  later  calculations. 

2.  The  Partition  Function  to  Second  Order 

The  trace  of  Eq.  (lOO)  with  respect  to  the  states 
of  Eq,  (101)  Is  indicated  by 


N 


Tr[exp(-^H)]  =1^  2iT^®^P  ^7^^ 

(102) 


n=0  n^+n2=n 


m=0 


«0 
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We  now  propose  to  expand  Eq.  (102)  explicitly  In 
terms  of  the  operators  D  and  U.  Since  U  raises  the  num¬ 
ber  of  excitations,  n,  and  D  lowers  n,  and  because  the 
states  of  Eq.  (lOl)  are  orthogonal  in  n,  it  follows  that 
all  non-vanishing  elements  of  the  trace  must  have  the  num¬ 
ber  of  D  operators  equal  to  the  number  of  U  operators. 

Thus  for  zeroth  order,  we  have  no  D  or  U  operators.  The 
first  order  correction  vanishes,  and  the  second  order 
contains  one  D  and  one  U  operator.  The  third  order  van¬ 
ishes,  and  the  fourth  order  contains  two  D  and  two  U  op¬ 
erators,  etc.  Thus 


Tr[exp(-/^H)  ] 


X  „  J  <n,njn2l5^^exp(-/E) 
n=0  n, +no=n 

^  ^  (103) 


=  -  “o 


m=0  ®  ^  ^  ^0  p=0  ®  ^  ^  ^ 


in  which  the  zeroth  and  second  orders  have  been  indicated 
explicitly. 


The  solution  to  Eq.  (1G3)  in  zeroth  order  Is  identi¬ 
cal  with  that  of  the  Islng  model,  H  =  +1.  The  fact 

that  we  do  not  know  how  to  write  down  this  solution  in 
three  dimensions  will  not  present  a  great  deal  of  diffi¬ 
culty  in  the  low  temperature  expansion  to  he  derived  from 
Eq.  (103). 

In  order  to  develop  a  low  temperature  expansion  of 
Eq.  (103),  we  investigate  it  as  a  function  of  exp  (-^A/2), 
since  with  1/kT,  A>  0, 

11m  exp(-^^2)  ->  G.  (1G4) 

T^O  ' 


Since  the  states  of  Eq.  (IGl)  are  also  eigenstates  of 
Hq  +  I,  we  may  re-express  the  zeroth  order  of  Eq,  (103) « 


Tr[exp(-^H)]^^^^^j^  order 


(105) 


n  +n  I  2  TT 1  f  E?^-I^ ^  ^  dE|n,n^n2>, 


n=0  n^+n2=n 


where  E^^  =  E^  +  n^y2.  Evaluation  of  the  Integral  yields 

order  ' 

+  exp(-|^A/2)  I  <l,n^n2|exp(-^I)  |l,n,  n2> 

ni+n2=l 


(106) 
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+  exp(-^A)^^^^X^^^2,n^n2|exp(-^I)|  2,n,n2>  +  ,..j, 

where  the  expansion  In  powers  of  exp(-^($l<V2)  has  been  made 
explicit. 

In  the  second  order  of  Eq,  (103) ^  we  may  also  re¬ 
express  the  terms  In  a  form  similar  to  Eq.  (105),  but  now, 
for  convenience  of  expression,  it  Is  desirable  to  place 
subscripts  on  the  operator  I  which  denote  the  sequence  In 
which  all  operators,  D,  I,  and  U,  are  applied.  That  is, 

Iq  means  that  I  is  applied  to  the  original  eigenfunction, 
|n,nj^n2>.  Is  applied  after  the  operator  U  is  attached 
to  jn,nj^n2>.  !_]_  is  applied  after  the  operator  D  is  at¬ 
tached  to  jn,nj^n2>,  etc.  Thus  we  have  for  the  second  or¬ 
der  correction. 


i  I  <n,n,n2|2^^exp(-/^E)dE  (107) 

n»0  n^+n2=n  t  j.  / 


n, n^n2> 
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+  exp(y3^2)  I  <l,nj^n2j[- 


exp(-^Q) 


(A+I^-Iq)^  exp (  ~^Iq)  3 1  l,n^n2>  +  . . .  , 


n=0  nj^+n2=n 


(108) 


^  -  ^0 


-  E;j^_2  -  ^  -  Iq  ^ ^1^2^ 


=  exp(-^o)/  I  <2,n,n2|u  exp ( -/^  ) ]>^-— i— 

(nj^+n2*2  (A+Iq-Ii) 

|2,n,n2> 

+  exp(^^2)  I  <3,n^n2|U  exp(-/^I_^)D - ^ - 2 


1  2 


3,n,n^>  +  . .,  , 


where  we  have  carried  out  the  expansion  to  the  first  power 
In  exp(-/^A/2) . 

Thus,  collection  all  terms  of  Eqs.  (106),  (107),  and 
(108),  and  using  the  relationships 


I  0>  =  0;  I  |l,nj^ng>  =  0, 


(109) 


which  result  from  the  fact  I  measures  the  Interaction 


85 


between  Individual  spin  deviations  from  the  Neel  state^, 
and  consequently  requires  at  least  two  deviations  to  give 
a  finite  measure.  Eq.  (109)  will  allow  us  to  drop  the 
subscripts  on  I  In  the  expression  for  the  partition  func¬ 
tion  belowj  since  there  remains  only  one  position  In  which 
I  gives  a  non-zero  result.  Thus  we  have  for  the  partition 
function  In  the  low  temperature  limit,  and  to  second  order 
In  D  and  U, 


Tr[exp(-/^H)  ] 


to  second  order 


(no) 


=  exp 


(-/Eo)|l  . 


+  1  <2,n,n5|uD — - — ^\2,n,n^>  +  ...\ 

ni+n2=2  2  j 

+  exp[-^(E^4^)  iJn  +  Z  <l,n,np}[-  D - pU 

'  ^  I  nT+no=n  A+l)'^ 


1  2 


'^^(Z^+I)^^l^"”l”2^  n  +n  =3^^" ^1^2  1^^  (A+I)^  3,n^n2> 

+  ...| 

+  exp[ -^(Eq+^  ]i  Z  <2,n,n2  |exp(-/6’l)  I  2,nj^n2>  +  ...1 
Pl+n2'^2  j 


•  •  • 


Eq.  (110)  may  be  simplified  by  noticing  thats 


<o|,D — 3^-ul.0>  =  I  <2,n,no|uD - |2,n,np>,  (ill) 

'  /  A  I  T  \  *i  ^  i  ^  *  t ^  I  -r  \  ^ 


(A+I) 


ni+n2''2 


and 


Z  <l5n„n„iD  — ■^=~— U  jl,nTno> 

m,=l  ■  ^  (a+I)""  ^  ^ 


il+n2 


Using  Eq,  (ill),  we  get  for  Eq,  (110 ), 


Tr[exp('^H)  second  order  (112) 

=  exp(-/5EQ)|l  +/5<o|.D  ■^^u|o>  + 

+  exp[7?(EQ+f)]|N  TOT’ 

|l,n,n2>  +  . .  .| 

+  exp  [  (  Eq+a)  ^  ^<2,  n^ng  [  exp  ( -/l )  |  2,  n^ng? 

I 


rU 


1  2 


+  .  •  • 


+  ... 


3.  Relationship  Between  the  Partition  Function 
Expansion  and  the  Ground  State  Energy  Series 

Examining  the  first  term  In  Eq,  (112),  we  note 
that  It  appears  to  be  a  power  series  expansion  of  an  ex¬ 
ponential  In  which  the  argument  Is  Just  the  ground  state 
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energy  perturbation  series  derived  In  Chapter  IIIj,  Eq. 

(32). 

<0|  exp[-^(EQ-D^-|;^yU  +  o..)]|o>  (113) 

=  <0|  exp(-/SEQ)[l+/^[>^^^U  +  ..,](0>, 

In  Appendix  we  examine  the  fourth  order  perturbation 
terms  of  Eq,  (102)^  and  find  that  Eq.  (113)  Is  valid  to 
this  order  also. 

ThuSj  assiimlng  the  complete  validity  of  Eq.  (113)  to 
all  orders^  the  term  which  In  zeroth  order  corresponds  to 
(-/SEq),  goes  over  the  exp ( state^' 
be  expected, 

B,  The  Excited  States 

1,  Perturbation  Treatment  for  the  General  Case 

Examining  the  terms  In  Eq.  (112)  which  In  zeroth 
order  correspond  to  exp[(y3(EQ  +  we  find 

^  2  "  ®(A+l)^  ^  jl,nj^n2> 

(114) 

*^l™2“-‘' 

+  . . .  \|l,nj^n2>. 


under  assumptions  similar  to  those  leading  to  Eq,  (113) 


In  Appendix  G,  Eq.  (ll4)  Is  also  shown  to  be  valid  through 
fourth  order, 

Eq,  (ll4)  does  not  make  plain  the  relationship  of  the 
excited  states  to  the  ground  state,  but  an  explicit  calcu¬ 
lation  will  serve  to  do  so.  We  take  for  the  normalized 
eigenstates  |l,n^n2>, 

|l,l,0> 

ll,0,l> 


Substitution  of  Eq.  (115)  into  Eq.  (Il4)  yields 

I  <l,n3_n2|exp[-^(EQ  +|)][l  +  +  . .  .]|  l.n^L^'s 

n, +no=l  l  y 

^  ^  (116) 

=  axp[-^(Eo  +f)l|N  +  4  ^ 


where  we  have  made  use  of  the  Kronecker  delta  relation¬ 
ships. 


Iexp(i(A-/^').j]  = 

^  exp[l(A-A' )*4:]  = 


(117) 


The  validity  of  the  binomial  expansion  of  the  term  involV“ 
ing  the  operator  I  depends^,  of  course^,  upon  the  condition 
that  the  expectation  value  of  I/a  have  a  magnitude  of  less 
than  unity.  For  the  states  involved  in  Eq.  (114)^  this 
condition  will  be  fulfilled  except  for  the  linear  chain 
with  spin  one -half,  which  we  treat  separately  later. 

The  symmetry  of  the  operators  D,  I,  and  U  with  re¬ 
spect  to  the  sublattices  indicate  that  both  terms  of  the 
sum  over  the  reciprocal  lattice  in  Eq.  (ll6)  are  equal. 
Therefore,  we  calculate  only  one  of  these.  For  the  effect 
of  the  operator  I  on  USA^:jo>,  we  have  by  means  of  the  def¬ 
inition  of  I  given  in  Eq.  (20), 


(-I)”'sj[^|G<  =  (2J)”'US"  1,0> 


(118) 


+  (2J)"‘(2”-l)J(l-a)(f)  I 


Insertion  of  Eq.  (ll8)  into  Eq.  (ll6),  with  the  subse- 
quent  commutation  of  operator  through  all  operators 

until  it  operates  on  and  the  final  evaluation  of  the 

resulting  terms,  yields 


I  _  <l,n^n2|exp[-^(EQ  +  f )  ]|l  +  ...||l,n 

^  ^  (119) 


l’^2^ 


Reversing  the  process  indicated  in  Eq,  (ll4)  allows  ns  to 
write  for  the  energy  of  the  first  excited  states^ 


■'first  excited  states 


<DU> 


(120) 


+  ^[l  + 


2(2S)'q-1; 


2(S^0-l)  A  V  - 


Eq.  (120)  indicates  that  there  are  N  first  excited  states^ 
each  of  which  lies  near  the  ground  state,  as  indicated  by 
the  presence  of  the  ground  state  energy  perturbation  se¬ 
ries  appearing  in  Eq.  (120).  The  Incremental  difference 
between  the  ground  state  and  the  first  excited  states  is 
given  by  the  terms  proportional  to  ^2  in  Eq.  (120). 

There  are  N/2  different  Increments,  each  characterized  by 
a  different  A.  Thus  each  X  corresponds  to  a  doubly  degen¬ 
erate  energy  level.  The  dispersion  relation  for  the  first 
excited  states  in  this  second  order  approximation  is  then 
given  by  the  relationship  between  the  Incremental  energy 
and  the  reciprocal  lattice  wave  vectors  A« 


(121) 


AE. 


'first  excited  states 


2^^  *  2(22 Kq-I)  ■  “2(S!(o-i 


+ 


2.  Comparison  with  Spin  Wave  Theory 

The  first  excited  state  energy  derived  by  means 

11 

of  the  spin  wave  theory  is  given  by 


AE 


'first  excited  states 


Examination  of  Eq.  (121)  for  large  splnp(the  condition  of 
validity  for  the  spin  wave  theory)^  reveals  Eq.  (121)  to 
be  the  first  two  terns  of  a  binomial  expansion  of  Eq. 
(122).  This  is  precisely  the  kind  of  relation  which  we 
found  for  the  comparison  of  the  perturbation  and  spin 
wave  treatments  of  the  ground  state.  There  is  an  impor¬ 
tant  difference^  however^  which  affects  the  usefulness 
of  Eq.  (121). 

The  relative  merits  of  the  spin  wave  approximation 
versus  the  perturbation  technique  for  the  determination 
of  the  low  temperature  properties  of  real  antlferromagnets 
depend  upon  the  amount  of  anisotropy  present.  The  very 
nature  of  the  perturbation  technique  indicates  its  appli¬ 
cability  to  antlferromagnets  with  strong  anisotropy.  The 
zero  order  Hamiltonian  represents  the  complete  suppression 


of  any  contribution  from  the  off-sixls  components  of  the 
spln^  and  In  order  to  represent  the  isotropic  Hamiltonian^ 
the  contribution  from  the  off-axis  components  must  be  com¬ 
pletely  restored.  In  the  ground  state  problem.  It  was  un¬ 
necessary  to  go  this  far  In  order  to  achieve  an  excellent 
representation  for  the  energy,  and  relatively  good  repre¬ 
sentations  of  the  short  and  long  range  order.  In  the  case 
of  the  excited  states,  however,  because  of  the  Boltzmann 
factor,  greater  importance  Is  attached  to  the  lowest  en¬ 
ergy  eigenstates  of  Eq.  (121),  and  It  Is  just  these  states 
which  are  most  affected  by  the  presence  of  an  anisotropy. 

A  finite  anisotropy,  no  matter  how  small,  removes  the  de¬ 
generacy  between  the  lowest  of  the  excited  state  energies 

and  the  energy  of  the  ground  state,  thereby  Introducing  an 
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energy  gap.-"^ 

The  ground  state  Is  adequately  represented  by  a  few 
orders  of  the  perturbation  because  Its  deviation  from  the 
Neel  state  Is  determined  by  a  sum  of  unweighted  contribu¬ 
tions  from  all  of  the  modes  of  the  lattice.  The  term  "un¬ 
weighted"  is  used  here  to  Indicate  the  absence  of  a  Boltz¬ 
mann  factor  In  the  ground  state  sum.  The  relatively  poor 
representation  of  the  adjustment  In  the  ground  state  ener¬ 
gy  via  zero-point  fluctuations  associated  with  a  long 
wavelength,  which  give  only  a  small  contribution  to  the 
total  adjustment  of  the  ground  state  energy  because  of  the 
nearly  perfect  alignment  of  spins  In  a  long  wavelength 


mode,  are  easily  compensated  for  by  the  much  more  adequate 
representation  of  the  short  wavelength,  high  energy  zero 
point  fluctuations.  This  property  of  the  perturbation  ex¬ 
pansion  is  of  little  avail  when  it  comes  to  determining 
the  low  temperature  behavior  of  an  antiferromagnet,  where 
it  is  the  low  energy,  long  wavelength  modes  which  are  most 
Important, 

The  Important  properties  of  both  the  perturbation  and 

spin  wave  treatments  of  the  excited  states  can  best  be 

seen  when  compared  v/ith  the  one  exact  treatment  for  these 
iin 

states.  This  may  be  done  as  soon  as  our  treatment  of 
the  general  case  is  extended  to  the  case  of  the  linear 
chain  with  spin  one-half. 

3.  Perturbation  Treatment  for  the  Linear  Chain 
with  ^In  (Me-Half 

We  now  return  to  the  development  of  our  perturba¬ 
tion  series  at  the  point  of  Eqs.  (107)  and  (108),  Here 
we  must  regard  the  fact  that  for  some  states  of  the  linear 
chain  coupled  with  spin  one-half,  all  three  energy  de¬ 
nominators  are  equal.  This  is  reflected  in  the  fact  that 
the  binomial  expansion  used  in  Eq.  (ll6)  is  not  valid  for 
these  states,  and  consequently  the  expression  in  Eq.  (119) 
cannot  be  evaluated  for  S  =  1/2,  Xq  =  2.  Instead  of  tac¬ 
itly  assuming  that  the  middle  energy  denominator  of  Eqs, 
(107)  and  (108)  differs  from  the  other  two,  we  evaluate 
the  Integrals  in  the  following  manner. 


n=0  nj^+n2=n 


(123) 


E  -  ^  -  A  -  ^  -  I^E  -  ^ 


N 


-^^ex;^  (75E)dE 


X 


n^  21T1. 

^  E  -  [Eq  +  .(AtllA  +  i^]^E  - 


N  f  exp[-/g(EQ-l-^-|^^4I 


exp[-/g(EQ-l-^-|^^4I^)] 

2  ”^^0  ^  ^ 


-exp[-/g(EQ+^-gQ)] _ 

[I^  +  -is^  +  -  (Ji\5-(lf +Io )  1  ^ 


[4,  + 


-/exp  [  -/i  (  Eq+^^+Iq  )  ] 

(E  +.(?A)<^ 
0  2 


nA 

2 


+  I 


+I 


l)l 


j|n,n^n. 


which.  In  the  limit  as  +  A  Iq>  i.e.,  when  all  three 
energy  denominators  are  equal. 


N 

=  11m  £ 

Iq-*Ij^+A  n=0 


r,  j-TT  _„<n,n,no 
ni+n2=n  •’12 
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njnj^n2>. 


In  place  of  Eq.  (108)^  we  get 
N 

11m  T  y  y  <n,n,np|UD  (124) 

Iq+A*"^-!  n=0  nj^+n2=ai 


X 


rij  ® 


The  term  with  n  =  1  In  Eq.  (123),  and  the  term  with 
n  =  3  In  Eq.  (124)  provide  us  with  the  entire  contribution 
which  Is  proportional  to  exp[ -^(Eq+^2)  ] .  In  the  limits 
Indicated  In  Eqs.  (123)  and  (124),  this  contribution  van¬ 
ishes.  The  vanishing  of  this  term  corresponds  to  the 
omission  of  the  term  In  Eq.  (119)  which  Is  proportional 
to  (A-4J)~^.  Thus  we  arrive  at  the  following  Instead  of 
Eq.  (121) 

AE  =  2J[1  +  (1-a)^  -  (l-a)^_?A  ].  (125) 

TT 

In  Figure  2,  we  show  the  dispersion  curves  for  the  spin 
one-half  coupled  linear  chain  In  this  approximation,  the 
spin  wave  approximation,  and  the  exact  solution.  From 
Figure  2,  It  Is  apparent  that  Eq.  (125)  best  approximates 


the  exact  solution  in  the  short  wavelength  region^,  l.e,^ 
A|6 1  =  rr/2.  This  bears  out  what  was  previously  said 
with  regard  to  the  treatment  of  the  ground  state.  For  the 
treatment  of  the  excited  states^  however,  the  large  gap 
which  Is  evident  in  Figure  2  between  the  lowest  states  and 
the  ground  state  in  this  order  of  the  perturbation,  pre¬ 
cludes  the  use  of  this  theory  to  obtain  meaningful  low 
temperature  thermodynamic  behavior  for  real  antlferro- 
magnets.  Higher  order  corrections  are  not  likely  to  close 
this  gap  rapidly,  in  the  sense  that  the  presence  of  even 
a  small  gap  profoundly  modifies  the  low  temperature  be¬ 
havior. 


V.  SUMMARY 


Perturbation  series  for  the  ground  state  parameters 
of  the  Heisenberg  antiferromagnet,  i.e.^  the  energy^  short 
range  order^  and  long  range  order^  have  been  presented  in 
a  quite  general  form  through  fourth  order  in  the  perturba¬ 
tion  Hamiltonian,  The  final  zero  order  Hamiltonian  is  the 
Ising  model j  H  =  final  perturbation 

Hamiltonian  consists  of  the  off-axis  components  of  the  in¬ 
teraction.  The  problem  of  the  N-body  divergence,  present 
in  any  straightforward  Raylelgh-Schrodlnger  perturbation 
treatment  based  on  the  Ising  model,  has  been  handled  in  a 
unique  way.  The  choice  of  an  initial  zero  order  Hamil¬ 
tonian  other  than  the  Ising  model  makes  possible  the  gen¬ 
eration  of  ordinary  Raylelgh-Schrodlnger  perturbation 
series  in  which  the  N-body  divergence  is  cancelled  in  each 
order,  thus  satisfying  the  requirement  that  a  systematic 
method  be  used  to  eliminate  non-physical  contributions 
from  the  series.  After  the  N-body  divergence  problem  has 
been  handled  using  the  initial  zero  order  Hamiltonian,  a 
way  is  found  to  predict  certain  infinite  classes  of  terms 
belonging  to  the  original  perturbation  series.  The  inclu¬ 
sion  of  these  terras  in  the  perturbation  correction  is 
shown  to  be  equivalent  to  a  shifting  of  the  definition  of 
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the  zero  order  Hamiltonian  to  the  Islng  model o 

The  perturbation  corrections  to  the  energy short 
I’ange  order^  and  long  range  order  of  the  Neel  state  are 
exhibited  as  explicit  functions  of  the  spin  and  number 
of  nearest  neighbors j,  and  also  of  an  anisotropy  parame¬ 
ter  which  essentially  gives  the  order  of  the  perturbation 
series.  For  comparison^  equivalent  expressions  are  de¬ 
rived  for  the  standard  spin  wave  treatment  of  the  ground 
state.  This  procedure  allows  a  complete  analysis  of  the 
relationship  between  the  two  methods.  On  the  basis  of 
this  analysis,  it  is  found  that  the  present  method  is 
superior  to  the  spin  wave  method  for  the  treatment  of  the 
ground  state  short  range  and  long  range  order  parameters 
in  all  one  dimensional  problems,  and  that  it  holds  an  ad¬ 
vantage  also  in  the  treatment  of  the  spin  one -half  cases 
in  two  and  three  dimensions. 

Experimental  determinations  of  the  sublattice  mag¬ 
netization  in  real  antiferromagnets  show  a  slightly  larger 
value  for  the  sublattice  magnetization  of  the  ground  state 
than  that  predicted  by  either  the  perturbation  technique 
presented  in  this  dissertation,  or  the  spin  wave  theory. 
Though  the  results  of  the  perturbation  treatment  are  clos¬ 
er  to  the  experimental  findings  than  are  the  results  of 
the  spin  wave  theory,  complete  confidence  cannot  be  placed 
in  the  assumption  that  higher  orders  of  the  perturbation 
will  not  further  lower  the  perturbation  results.  On  the 


basis  of  the  analysis  of  the  relationship  between  the  two 
methods^  it  is  reasonable  to  believe  that  the  exact  sub- 
lattice  magnetization  for  the  Heisenberg  model  ground 
state  lies  somewhere  between  the  spin  wave  prediction  and 
the  perturbation  prediction. 

The  treatment  of  the  first  excited  states  of  the 
Heisenberg  antiferromagnet  by  means  of  the  perturbation 
theory  show  the  anisotropy  to  have  a  strong  Influence  up¬ 
on  the  energies  of  the  lowest  of  these  states,  Thus^  be¬ 
ginning  with  a  zero  order  Hamiltonian  which  is  very  anlso 
troplcj  and  calculating  only  a  few  orders  of  the  perturba 
tlon  series^  an  adequate  representation  for  the  energies 
of  these  states  is  not  found.  It  is  concluded  that  one 
should  rely  on  a  continuation  of  the  spin  wave  treatment 
of  the  states  in  the  manner  presented  by  Oguchi”^  for  any 
Improvement  in  the  representation  of  the  low  lying  ex¬ 


cited  states. 
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Linear  Chain 


Plane  Quadratic 


Simple  Cubic 


FIGURE  1.  Classical  Antiferromagnetic  Ground  States 


TABLE  I 


Qualitative  Results  of  Previous  Methods  Used  to  Compute  the 


Sublattice  Magnetization 

of  the 

Ground 

State 

Is  there  a 
Finite  Sublattice 
Magnetization? 

Linear 

Chain 

Plane 

Quadratic 

Simple 

Cubic 

S=l/2 

S>l/2 

S=l/2 

S>l/2 

S«l/2 

S>l/2 

Anderson 

No 

No 

Yes 

Yes 

Yes 

Yes 

Kubo 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 

Kasteleljn 

No 

Taketa  and 

Nakamura 

No 

No 

Marshall 

No 

No 

No 

Fisher 

Yes 

Yes 

Davis 

Yes 

Yes 

Yes 

Yes 

Yes 

Yes 

TABLE  II 


The  Perturbation  Correction 

to  the 

C 

Energy 

2 

ijS-t/ibic© 

S=l/2 

S=1 

S=3/2 

S=2 

Linear  Chains  This  Work 

Linked  Cluster 

Plane  Quadratics  This  Work 
Linked  Cluster 

Simple  Cubic s  This  Work 

Linked  Cluster 

IcOOOO 

1.0000 

0,3333 

0,3333 

0,2000 

0,2000 

0,3333 

0,3333 

0 , 1429 
0,1429 

0.0909 

0.0909 

c 

0,2000 

0,2000 

0,0909 

0,0909 

0,0588 

0,0588 

'4 

G„1429 

0,1429 

0,0667 

0,0667 

0,0435 

0,0435 

Linear  Chains  This  Work 

“0 , 2500 

0,0782 

0,0591 

0,0456 

Linked  Cluster 

-0,4590 

0,0292 

0,0284 

0,0226 

Plane  Quadratics  This  Work 

“O0O065 

0  e  0251 

0,0099 

0,0078 

Linked  Cluster 

-0.0098 

0,0123 

0,0098 

0,0078 

Simple  Cubic s  This  Work 

-O.GOO8 

0,0055 

0,0045 

0,0037 

Linked  Cluster 

-0.0015 

0,0054 

E/Eq( 

0,0045 

a=0)* 

0,0037 

Linear  Chains  This  Work 

1,7500 

1.4115 

1,2591 

1,1885 

Linked  Cluster  I 

1,5410 

1,3625 

1,2284 

I0I655 

Linked  Cluster  II 

1,7363 

1.3567 

1,2287 

1,1680 

Plane  Quadratics  This  Work 

1.3268 

1,1680 

1,1008 

1,0745 

Linked  Cluster  I 

1,3235 

1,1552 

1 , 1007 

1,0745 

Linked  Cluster  II 

1.3281 

1,1563 

1,1027 

1,0765 

Simple  Cubic s  This  Work 

1,1992 

1,0964 

1,0633 

1,0472 

Linked  Cluster  I 

1,1985 

1,0964 

1,0633 

1,0472 

Linked  Cluster  II 

1,1996 

1,0973 

l,0643 

l,048l 

* 

The  greater  the  magnitude  of  this  term^  the  lower  the  cor¬ 
responding  energy since  Ej,  E^  <  0, 


104 


TABLE  III 

Comparison  of  the  Perturbation  and  Spin  Wave 
Corrections  to  the  Energy 


E/EQ(a=0) 


Lattice 

S=l/2 

S=1 

.  s=3/2 

S=2 

Linear  Chains  This  Work 

Corrected  Spin  Wave* 
Hybrid  Method 

1.7500 

1.8589 

1.7651 

1.4115 

1.3964 

1.4485 

1.2591 

1.-.2569 

1 0  2903 

1 . 1885 
1,1900 
1.2144 

Plane  Quadratics  This  Work 
Corrected  Spin  Wave 
Hybrid  Method 

1.327 

1.341 

1.322 

lc.l68 

1.164 

1.177 

1 . 101 

1.108 

1.109 

1.075 

1.081 

1.081 

Simple  Cubic s  This  Work 

Corrected  Spin  Wave 
E^brid  Method 

1.199 

1.203 

1.199 

1.096 

1.099 

1.099 

1.063 

1.066 

1,066 

1.047 

1.050 

1.049 

Spin  wave  approximation  +  first  order  correctiono 


TABLE  IV 


Comparison  of  the  Perturbation  and 
Spin  Wave  Corrections  to  the  Short  Range  Order 


_ -^(a=0) _ 

La.ttlce 

S=l/2 

S=1 

S=3/2 

S=2 

Linear 

Chain;  This  Work 

0.7500 

0.4321 

0.6227 

0.7203 

Plane  Quadratic;  This  Work 
Corrected  Spin  Wave 
Hybrid  Method 

0 . 686 
0.970 

1.038 

0.782 

0.796 

0.763 

0.879 

0.822 

0,822 

G.9I0 

0.850 

0.851 

Simple 

Cubic;  This  Work 
Corrected  Spin  Wave 
Hybrid  Method 

0.802 

0.856 

0.877 

0.907 

0.886 

0.903 

0.928 

00  915 
0.916 

0.945 

0.933 

0.933 

TABLE  V 


The  Second  and  Fourth  Order  Corrections  to  the 
Long  Range  Order 


Lattice 

S=l/2 

S=1 

8=3/2 

S=2 

_s^ 

Linear  Chain 

1.0000 

0 , 2222 

0.1200 

0.0816 

Plane  Quadratic 

0.2222 

0.0816 

0.0496 

0.0356 

0.0178 

Simple  Cubic 

0.1200 

0,0496 

0.0311 

D4 

0.0227 

Linear  Chain 

0.2500 

0.1433 

0.0687 

0.0458 

Plane  Quadratic 

0.0356 

0.0270 

0,0181 

0.0136 

Simple  Cubic 

0.0080 

0.0107 

0,0083 

p(a=0) 

0.0061 

0,0050 

Linear  Chain  -O.25OG  0,63^5  O.8II3  0,8726 

Plane  Quadratic  0.7422  0,8914  0.9323  0.9509 
Simple  Cubic  0.8720  0.9397  O.9606  0.9712  0,9772 


TABLE  VI 


Perturbation  and  Spin  Wave  Corrections  to  the  Long  Range 
Order  of  the  Isotropic  Ground  State 


Lattice 

S=l/2 

S=1 

s=3/2 

S=2 

S=5/2 

Linear  Chains 

This  Work 

-0.2500 

0.6345 

0.8113 

0.8726 

Plane  Quadratics 

This  Work 

0.7422 

0,8914 

0.9323 

0.9509 

CorrectedlSpln  Wave 

0.607 

0,803 

0.869 

0.902 

Hybrid  Method 

G.597 

0,801 

0.868 

0.902 

Overlap 

0.752 

0.893 

0.933 

0.951 

Simple  Cubic s 

This  Work 

0.8720 

0.9397 

0.9606 

0.9712 

0.9772 

Corrected  Spin  Wave 

0.844 

0.922 

0.948 

0.961 

0,969 

Hybrid  Method 

0.848 

0.922 

0.948 

0.960 

0.969 

Overlap 

0.868 

0.940 

0.961 

0.972 

0.977 

TABLE  VII 


The  Square  of  the  Long  Range  Order  Parameter  for  Comparison 
with  the  Short  Range  Order  Parameter 


Lattice 

S=l/2 

- r— 

S=1 

S=3/2 

S=2 

Linear  Chain;  This  Work 

0.403 

0.658 

0.761 

Plane  Quadratic;  This  Work 

0.551 

0.794 

0,869 

0.904 

Simple  Cubic;  This  Work 

0.760 

0.884 

0,924 

0.943 
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APPENDIX  A 


In  this  Appendix  we  calculate  the  sixth  order  energy- 
correction,  For  the  general  case  with  arbitrary  spin  and 
lattice^  we  proceed  as  far  as  the  sixth  order  equivalent 
of  Eq,  (52)5  Chapter  III.  For  the  case  of  the  linear 
chain  with  spin  one-half^  we  continue  on  to  determine  ex¬ 
plicitly  the  sixth  order  energy  correction^  which  is  found 

to  be  identical  to  that  found  for  this  case  by  a  different 
IQ 

method.  ^  The  desirability  of  such  an  effort  lies  with 
its  implication  that  we  have  indeed  recovered  the  Islng 
model  as  a  zero  order  Hamiltonian^  at  least  through  six 
orders  of  the  perturbation. 

We  begin  by  considering  Eq.  (53) >  Chapter  III. 
applying  the  relationships  indicated  in  Eq.  (4o)  and  (51)^ 
the  operators  I  which  are  superscripted  with  the  symbols 
P;,  r,  s,  and  q  may  be  replaced  by  the  explicit  expression 
of  their  effects  upon  the  U  and  D  operators  of  the  matrix 
elements.  In  performing  this  replacement,  we  are  aided 
by  the  Hermltlan  property  of  the  matrix  elements,  in  par¬ 
ticular  by  the  fact  that  <o|diPdi^1  =  (|  I^UI%|o>)'^,  Fol¬ 
lowing  this  replacement,  the  sums  over  p,  r,  s,  and  q  may 
be  performed  to  give  in  place  of  Eq.  (53)* 
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Here  we  have  subscripted  the  D  operators  with  the  sub¬ 
scripts  77'  and  88'.  These  have  an  entirely  equivalent 
meaning  with  respect  to  the  D  operators  as  the  11 '  and 


22*  subscripts  have  with  respect  to  the  U  operators.  Sub¬ 
scripts  utilizing  the  intermediate  numbers  33'  to  66*  will 
be  needed  when  the  operator  is  replaced. 

In  order  to  perform  the  sum  over  m^  we  need  to  know 
the  effect  of  l”’*  on  U^|o>j  UUj^U^,|o>j  inJ^2,U2^i2|0>,  and 
DU^|o>,  ®%2*^l'2l*^^*  detemlned  in  a 

manner  analogous  to  that  used  to  determine  the  effect  of 
the  operator  on  U^jo>,  that  is,  by  applying  successive¬ 
ly  higher  powers  of  the  operator  I  to  each  term.  As  we 
saw  in  the  series  of  equations  leading  to  Eq.  (51)^  i.e., 
Eqs.  (42)  and  (44),  each  application  of  the  operator  I 
serves  to  Introduce  new  terms  derived  from  the  old  by  link¬ 
ing  two  operators  with  a  new  subscript,  as  in  Eq.  (42)  and 
Eq.  (44),  When  all  operators  have  been  linked  with  a  sub¬ 
script,  e.g.,  U^2i%i2l‘^^^  (^5)^  the  resulting  com¬ 

pletely  linked  function  is  an  eigenfunction  of  I,  and  the 
production  of  new  terms  ceases  upon  further  application 
of  I.  There  are  sixteen  different  operators  which  can  be 
formed  from  the  basic  operator,  U^,  by  the  linking  process. 
We  list  them,  utilizing  an  extension  of  the  subscript  no¬ 
tation  begun  when  the  different  linked  operators  formed 
2 

from  U  were  defined.  For  clarity,  we  follow  the  list  with 
the  explicit  expression  for  some  of  these  subscripted  op¬ 
erators,  the  definition  of  the  others  following  directly. 


(A2a) 

u3. 

(A2b) 

UUiUj^,, 

(A2c) 

UU^2-^1-2^ 

(A2d) 

^2^1 ^1 ' 2 ' ^ 

(A2e) 

(A2f) 

U2Ui2,U,„ 

(A2g) 

^3^12 '3'% '2" 

(A2h) 

U3 ,  U^2 ' 3%  '  2^ 

(A21) 

^23'%3^1'2'" 

(A2J) 

^23 '^12'%  >3" 

(A2k) 

^34'%2'3'%'24" 

(A21) 

^34^12'3'^1'24*^ 

(A2m) 

^34'%2'3'4^1'2^ 

(A2n) 

^3'4'^12'3%'24" 

(A2o) 

^34'5'^12'3'5%'24" 

(A2p) 

^34'5'6^12'3'5’^1'246 

The  definitions  of  lIUj^U^,  and  being  simple  exten¬ 
sions  of  U^Uj,  and  U^2'%'2'’  definitions  were  given 

In  Eqs.  (43)  and  (45),  respectively,  we  begin  with  the 
definition  of  1 2 1  • 


(A3a) 
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==  (\)  (\)  (s-ir^-fg)!  °>' 

! 

U^Uij.UiJo  -  (|)2j3(i.a)3^^I^^  V^'^iW  («o) 

X  (s!^) (s-^)  (^,^.5^)1  0>, 


'’3“l2.3.''l.2l<»  ■  («<») 


^  (X+^)  (2^2>  (2/ij+Si)  ‘^2-^23^^  <\t€i>  <^-<12-63) 


"34'5'6''i2'3'5”i'246'I°>  = 

(A3e) 

^ 

We  are  now  In  a  position  to  express  the  application 
of  l”*  to  U^l  0>,  UU^U-j^,|o>,  and  TJU^gt^l  1 2  terms  of  a 

linear  combination  of  functions  made  up  from  the  operation 
of  Eqs.  (A2)  to  the  Neel  state,  Ey  successive  application 
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of  the  operator  I,  we  determine  that. 


=  (-2j)®|3®a3|o>  +  6(4"^-3®)uu^U3_,  |o>  (a4) 

+  3(5®-2.4”'+3“)(UU3^2,Ui,2|0>+U2U3^U3^,2,  |0> 

+U2 ,  .  I  0>+2U2U^2  •  ^1 '  I 

+  2  ( 6"'-3 .  s'^+s .  4“-3”')  (  3U3U^2 ' 3 '  ^1 ' 2I 

+3U3,Ui2,3U;^Jo>+3U23, 

%3%*2»  I  °^'^^23'%2'%'3l°^^ 

+  3(7”‘-4.6"'46.5™-4.4"'+3®)  x 

+  U32j.U22'3*%' 24' f  ^34 '^12' 3 '4^1 '2!®^ 

+  ^3'4'^12'3^1'24|°^^ 

+  e{Q^-5, 7^+10 . 6*"-10 . 5®+5 . 4'®-3“^) 

X  (U3^ 

'5'^12'3'5%'  241  o» 

+  ( ^-6 . 8"‘+l5 . 7®-20 .  e^+is .  5^-6 . 4“+3"^) 


^  ^^34'5'6^'12'3'5%'246' 


and  for  the  operation  of  l”*  onto  UTJ^U^jo,  we  determine 
that 
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+U2 ,  ,  1 0>+2\J^V^^ ,  ,  1 0>) 

+  (6’"-2 . 5"*+4™)  (  ' 3 '  ^1  •  2  |o>+3U3  ,  ,  3U^ ,  gl  0> 

+3^23  •  ^13^1  '  2 '  I  0>+U23  I  U]_2 '  ^1 ' 3 

+  2  (  t'^-S  .  6"'+3 . 5“-4"*)  ( 2U34 .  %  2 . 3 .  Ui .  24 1 0> 

+U321U12 1 3 1  1 24 1  I  ^^■^^34 1  ^X2 ' 3 ' 4^1  *  2 1 

+U3 , 2j^XI^2  1 3%  124!®^^ 

+  5  ( 8“-4 .  r'^+e .  6^-4 , 5”^-h4”‘)  U3^ ,  ^ ,  U3^2  •  3  •  5^1 . 24 !  °> 

+  (9"^-5.8®+io.7®-io.6’^+5.5'"-4’^) 

^  (^34'5'6%2'3'5^1'246' 1°^^} 

Finally,  for  the  operation  of  to  UUj^2,U^,g[o>,  we  get 

l"'uUj^2'^l'2i°>  "  (-2J)“|5’”uU;^2'^^'^^°^ 

+  2 ( 6^-3^) ( U3U^2 • 3^1 ' 2I  °>+^3 • ^12 ' 3^1 ' 2 1 

+  ( 7'"-2 .  e^'+s”)  ( 2U34 .  U12 . 3 .  ,  24  |o> 

+U34U12 1  2 • ' 24 ’ f  0>+2U3^ ' ^12 ' 3 ' 4^1 ' 2  f 
+^314^^12 '3^1 '24!°^^ 

+  4  ( s^^-s .  7"'+3 . 6^-5^ )  U3^ ,  5 ,  2' 3 ' 5^1 ' 24| 

+  (9“-4.8“+6.7®-4.6™+5"*)U34,3,6U^2'3'5^1'246'I  0>}- 
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For  the  successive  application  of  I  to  the  functions 
DU^|o>,  DU^Uj^ijo,  and  , 2  |o>,  we  cannot  use  the 

definition  of  I  given  In  Eq.  (20),  Chapter  III,  since  It 
was  defined  only  in  terms  of  Its  application  to  U-llke  op¬ 
erators.  If,  however,  we  first  commute  the  D  operator 
through  all  U-llke  operators  until  It  operates  on  the  Neel 
state,  we  are  left  with  functions  which  Involve  only  U-llke 
operators,  l.e.. 
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)(S 


)]. 


In  order  to  simplify  the  expression  of  the  results 
which  follow,  we  introduce  the  following  abbreviation. 


DU' 


=  [2<Q>u|o>  +  R'|0>]. 


(A9) 


For  successive  applications  of  I  to  1 1 0>,  we  get 

i“dU^U3^,|0>  =  I™[2(|)Kq<Q>u|o>  +  R^^,|o>]  (AlO) 

=  DU^Uj^,|o>j  m  =  0 
=  (-2J)"*DUjUj,|  m  0 


Since  all  U-like  operators  are  maximally  connected 
in  DUj^2i'^1!2  1®^^  ’”'® 

l'"l>Ui2,Ui,2|0>  =  (-2J)"'DU3^2'Ul'2f°>- 

The  insertion  of  Eqs.  (a4)--(A7),  and  Eqs.  (AlO)-(All) 
into  Eq.  (Al),  with  the  subsequent  sumiTiation  over  the  in¬ 
dex  m,  results  in  the  sixth  order  analog  of  the  fourth 
order  expression  given  in  Eq.  (52),  Chapter  III.  Thus  we 
have  completed  the  first  phase  of  the  Appendix. 

The  second  phase  consists  of  determining  the  explicit 
sixth  order  energy  correction  to  the  linear  chain  coupled 


spin  one-half.  The  first  step  in  this  operation  is  to 
write  out  the  results  of  the  summation  over  the  index  m 
for  thpse  matrix  elements  which  contribute  to  this  case. 

If j  after  the  summation  over  the  index  m,  a  matrix  element 
has  a  coefficient  Involving  the  factor  (l-4j/2i)  then 
the  corresponding  matrix  element  must  be  null  for  the  case 
of  the  linear  chain  with  spin  one -half.  If  this  were  not 
the  case,  then  the  perturbation  series  would  diverge^ 
since  (1-4j/A)  =  0  for  this  case.  We  prove  this  for  two 
matrix  elements  in  Appendix  B,  and  assume  it  for  the  rest. 
The  summation  over  the  index  m  for  the  terms  contributing 
to  the  sixth  order  energy  correction  under  the  above  as¬ 
sumption  thus  yields 


123 


8/3  (  J/a)  ^  '  2 '  <^7^7 '  ^^2 '  ^1 2^1 '  >R 

_ +  2<Dy Dy I ^^2^1 2 ' ^1 ’ ^R _ 


(i-ioj/3a)  (i-3J/a)  (i-8J/3a)  (i-2J/a) 

^  |<D^UDU%J^  2(J/a)<D^UDuSj^^^^ 
[(1-2J/a)^  ( 1 -2 J/a  ) 5 


4  (  J/a  )  <D^UDU3^Uj_  ,  >j^  8  (  J/a)  '  >R  ^  ^  ^ 

(1-3J/a)(1-2J/a)'^  (1-3J/a)(1-2J/a)5 

4  (  J/a)  ^<^7^7 1  ■  >R  8  (  J/a  )  ,  UDUiUi;>r  ^  ^  ^  1 

^  (1-3J/a}^(1-2J/a)^  ^  (1-3J/a)^(1-2J/a)^  ' 


There  are  apparently  fifteen  matrix  elements  to  com¬ 
pute  in  Eq.  (A12).  We  can  use  the  Hermltian  property  of 
the  matrix  elements  to  reduce  the  number  of  elements,  since 
they  are  not  linearly  independent.  Using  the  relationship 
<oli)^I  -  (lU^jo>)'’',  along  with  Eq.  (A4)  with  m  =  1  and  Eq. 
(A5)  with  m  =  1,  we  obtain 


<d2iUUj^U^  ,  >j^  =  (  -2 J)  j^3<D2uU3^U^  ,  >j^  +  6<D^D^ ,  DUU^U^  , 
=  (-2J)|4<d3uu^U3^,>j^  +  <D^UU3^2'^1'2^R 


+  2<D^U2Uj^,2,U^,>j^  ,. 
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Por  S  =  1/2,  ITq  =  2,  Eg.  (A13)  yields 

<lPu2%%  ,  2 1  >R  +  <D^U2  ,  2Uj_ ,  >pj  +  2<d3u2U^  i  2 '  %  ' 


=  6<DyD,jr ,  DUUj^U^  J  >J^  -  <D^IIU^U2^  ,  >J^, 


(A14) 


since  UU]_2|Uii2^R  ”  ®  this  case.  This  Is  one  of 
the  matrix  elements  which  we  show  explicitly  to  he  null 
In  Appendix  B. 

Using  the  same  Hermltlan  property  and  a  combination 
of  Eqs.  (a4)-(A6),  we  obtain  the  following  two  equalities. 


<d3iU2UiUi.2.>r  +  +  2<I)3lU2U^2 '^1  •  >P 


=  (-2J)|5[<d2u2UiU3^,2,>r  +  <I>\.Ui2Ui.>r 

+  2<d3u2U^2'Ui->R^ 

+  2  [  2 ' 3 '  ^1 ' 2^R  2<D^U2  i  U^  2 1 2^1 ' 2^R 

+  3'<I5^U23  '  ^13^1  ’  2 '  ^R  *^^^23  ’  ^1 2  •  ^1  *  3^R 
=  (“2J)^3[>1^U2^1^1' 2'^R  ^^^2 ' ^12^1 ' ^R 

+  2<d3u2U3^2'^1'^R^ 


(AI5) 


since  for  the  case  of  the  linear  chain  with  spin  one- 
half  we  have  assumed  the  following  null  relationship, 

2<D^U2TJ2^  2 1 2  •  •  2^R  2<D^U2  i  2 1 3^1 » 2^R  ^  ^ 

+  3<D^U23,U^3U^,2i>r  +  <^^23'^12'%'3^R  “ 

we  obtain  for  this  case, 

<D^D^,DU2U3^U3^,2,>j^  +  <D^D^,DU2,Uj^gU3^,>R  (AI7) 

+  2<D^D^,DU2U^2'%>>R  =  ^/3 

=  l/3t  <B^2Di.0i2.>r  + 

4.  a<D3n2D,2,u,,>„]. 

with  the  reatlonshlps  Eq.  (Al4)  and  Eq.  (AI7),  the 
number  of  Independent  matrix  elements  Is  reduced  to  nine 
for  Eq.  (A12). 

The  first  step  In  the  calculation  of  these  nine  ma¬ 
trix  elements  Is  to  commute  each  D-llke  operator  through 
all  U-llke  operators  until  It  operates  upon  the  Neel  state. 
The  results  are  then  given  In  terms  of  the  commutators  be¬ 
tween  D  and  U  as  defined  In  Chapter  III,  Section  A,  Part  2, 
It  should  be  remembered  that  the  subscript  R  on  the  matrix 
elements  Indicates  that  all  terms  proportional  to  powers 
of  N  greater  than  one  are  to  be  disregarded,  as  shown  In 


Chapter  III,  Section  A,  Part  2,  Thus  we  get  for  the  nine 
matrix  elements  of  interest, 

=  2<R'’“R>  +  <V>,  (Al8a) 

<D^UUj^Uj^,>  =  6(|)yQ<P><Q>  +  +  <Rj|^Rj^,>  + 

(Al8h) 

<V7'“^lUl'>R  = 

+  <R^,^,R2^2^,>  +  <Rj^y^,Ru^,>  +  <V2.j^iyyi>, 

<D^UDU^>j^  =  <R'^R>,  (Al8d) 

<D^UDU3^Uj^,>j^  =  2(|)^q<P><Q>  +  <R+Rj^^,>,  (Al8e) 

<D^D^,UDU^U3^,>j^  -  4{|)yQ<F><Q>  +  +  <R+^,R^^,> 

(Al8f) 

<D^UDU^>^^)  =  <R’^R“>,  (Al8g) 

=  2(|)yQ<F><Q>  +  <R'^R“>,  (Al8h) 

<D^D^,UDU^U3^,4^^  -  4(|)yo<QxP>  +  4(|)2^2<q^  +  <R+^,R'> 

(Al8i) 
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(19c) 


To  complete  the  evaluation  of  these  matrix  elements, 
the  commutators  of  Eqns.  (iSa-l)  are  replaced  by  their  ex¬ 
pressions  In  terms  of  spin  operators,  being  careful  to  ob¬ 
serve  the  linking  of  operators  Indicated  by  the  subscripts 
The  resulting  general  expressions  are  then. 


A  V  <- 

^0 

-  216(SYq)5 

-  54(3^^)  + 


x^l2(SyQ) 


6  ( 15  -45^0-^40) 


^0 


(A19a) 


V  2  ^  Y  2 


^1' 


<d3uu. 


1^1, >R  =  (^)2^  X  jecs^o)^  (A19b) 

*0  ^ 


f-  6^3+18/ +49 
^  6(siro)«[-2 - 0  0  ; 

.  2T  3+16Y  -7 

-  12(SYo)5[-0 - 0-— 0_] 

^0 


+ 


6(sro)‘*[^4o  -  54(s)ro)3  +  9(sio)2  , 
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/+  +  321^3  +  .  351^^,  +  £03 

I  °  - 


.  ny  3  +  101/  2  _  46^  -  77 
-  2(s^„)5  - 2 - 

^0 


<D%IIU^>j;  =  (!!^)26(3j^)2 


(Aiga) 


:{t(S>o) 


4  (15Yo"-45iro44o)  _  ,3 


-  t 


,  _  (4K,^+3IC,-3)  1 

'r°  -  -  8(S«o)  *  ij- 


<D=UBOini.>K  =  (^)2® 

0 


(A19e) 


[2(3^0 


.  ,6  -  63/0  ■>■  49 


^  p.  r_2  +  er^  -  6  h  3^0^  +  23^-2 

-  8(s)ro)5  -2_0 - ^  6(3y^)4  ,0  p - 


-  8  (siro)3  +  (sYj,)^  }, 


<DyDj,UDUj^U^,>^ 


(A19f) 


j  2  '  c;  ' 

0 


X  V  (SYq) 


6.  +  24^^  +  18^^  -  99^  H-  67 


y; 


p-  -I-  SiTf.  ~  3  ^  h  SVo^  +  3^0  "  3 

-  16(8^0)5  +  4(SY^)^  --^—7^ - 

^0  ^0 


-  8  (S^)3  +  (3)^)2  j, 
<D^UDU^>;^^^  =  •— (^)2^ 


(A19g) 


X 


(36(3ir„)S  -  48(afo)5 


I  °'  >0 


+  4(s>o) 


4 


8(sKo)3  +  (sir^)^j. 


<D^UDUj_U3_,>^ 


,(i>  .  £1^ 

"0 


(A19h) 


x^i2(3y„)6 

'0 


.  R  3^0^  +  6  y  A  +  2/0-2 

-  8(31^0)5  -S—^O -  ^  6(s^^)4  -0.-^.^ - 


-  8(sYg)^  +  (sy^)^l 
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(A191) 


X  W)fi  _V  ^  ^  fo'  -  -  9 

I  /q 


i6(s^q)5  Jo._1.^^P 
^0 


i  .  yo. .-  ^ 

^0 


8(s!ro)3  +  (SY^)^}. 


Putting  S  =  1/2  and  Kq  =  2  Into  Eqs.  (A19a-1)  and 
using  Eqs,  (Al4)  and  (A17)j  provides  numerical  values  ifor 
all  necessary  matrix  elements  In  the  calculation  of  the 
sixth  order  energy  correction  for  the  linear  chain  with 
spin  one-half.  Substitution  Into  Eq.  (A12)  yields 

(l-a)^  IQ 

^S=1/2j^q=2  ^  agreement  with  Walker.  ^ 


APPENDIX  B 


In  this  Appendix,  we  show  explicitly  that,  for  the 
linear  chain  with  spin  one-half, 

=  0,  (Bl) 

and 

<d\t)U^2,U^,2>j^  =  0  (B2) 

Eqs.  (Bl)  and  (B2)  represent  typical  terms  resulting 
from  Eq.  (53) ^  Chapter  III,  which  have  denominators  which 
vanish  for  this  case.  It  should  be  realized  that  the 
vanishing  of  the  denominators,  coupled  with  the  vanishing 
of  the  matrix  elements,  does  not  Indicate  an  Indeterminacy 
In  the  evaluation  of  the  corresponding  term  In  Eq,  (53). 

In  the  course  of  computing  succeeding  orders  of  the  per¬ 
turbation  correction,  matrix  elements  like  Eqs.  (Bl)  and 
(B2)  appear  with  finite  coefficients  early  In  the  series. 
It  Is  at  this  point  that  their  null  property  should  be 
taken  Into  account,  and  not  after  the  sum  of  all  possible 
coefficients  throughout  the  perturbation  series  has  been 
determined  for  a  given  matrix  element  under  the  assiimptlon 
that  It  Is  not  mill. 
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We  consider  first  the  form  of 


<D  %1'22'^ 


(B3) 


2 '  ^1 ' 2^  "I  I PPi>^ 


11'22' 


which  is  determined  by  means  of  the  commutators  of  D  and 
U  defined  in  Chapter  III,  Section  A,  Part  2.  The  sub¬ 
scripts  refer  to  the  linking  of  operators  discussed  in 
Appendix  A.  Substitution  of  the  commutators  with  their 
expressions  in  terms  of  spin  operators,  and  the  subsequent 
evaluation  of  these  expressions,  yields 


■'^^12'%'2^R  V  P  (  p 


^0^ 


(-2^)2^ 


(B4) 


xl2l6(S^) 


I 


^  -  24(S^)5  -£^-11 


'0 


+  6(s^) 


4  10 +  27^0  -  27 


54(s/o)3  +  gisy^fj 


which  for  S  =  1/2,  =  2,  gives 


<D  I  2^R  ® 


(B5) 
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p 

For  the  form  of  <D  UDUj^2,U^,2>j^,  we  find  upon  apply¬ 
ing  the  commutation  relations  between  D  and  U, 

<D^UPU3^2«^1'2^R  "  '  22 


which,  \ipon  replacement  of  the  commutators  with  their 
equivalents  In  terms  of  spin  operators  and  subsequent 
evaluation  thereof,  yields 


<d2udU^2'Ui.2>r  = 


Substitution  of  S  =  1/2,  =  2  Into  Eq,  (B7)  then  yields 


the  predicted  result. 


APPENDIX  C 


In  this  Appendix^,  we  will  compute  the  fourth  order 
perturbation  correction  to  the  partition  function  as  It 
is  expressed  In  Eq.  (103)^  Chapter  IV.  The  computation 
Is  not  carried  out  to  the  point  of  providing  a  nmerlcal 
correction^  expressed  as  a  function  of  spin  and  lattice, 
but  only  so  far  as  to  show  the  validity  of  relationships 
such  as  Eq.  (113)  and  Eq.  (ll4)  to  fourth  order.  The 
fourth  order  perturbation  correction  Is  given  by  the  non- 
vanlshlng  elements  of 


n®0  rij  4-n2=n 


X 


The  non-vanlshlng  elements  of  Eq.  (Cl)  are,  of  course, 
those  In  which  the  nximber  of  D  operators  Is  equal  to  the 
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number  of  U  operators.  These  terms  correspond  to  the  fol¬ 
lowing  ordering  of  the  D  and  U  operators, 

DUDU,  (C2) 

UDUD, 

DUUD, 

UDDU, 

DDUU, 

and  UUDD, 

Evaluation  of  the  corresponding  six  Integrals  obtained 
from  Eq.  (Cl),  and  an  expansion  of  each  in  powers  of  exp 
(y^A/2),  yields  the  following  result  for  the  siun  of  con¬ 
tributions  from  all  six  integrals,  plus  the  zeroth  and 
second  order  corrections  of  Eq.  (110). 


Tr[e35>(-/5H)]j^  (03) 

=  exp(-y6EQ)-|l  +  <o|o^|o>  +  T.  <2,n^n2[  0g  ]2,n^n2> 

I  ni +no=2 


1  2 


+  I  <4,n^n2(0cl4,n^n2>] 

^1™2~^ 

+  exp[-^(EQ+V2)[N  + 

4*  ^  I  1 

+  „  I  -=  <5,n^n2|0  |5,nin2>]  + 

n-,  +nr,-h 


1  2 


•  9  •  J 


where 


0^  =  [-  1>— -  oU  +  2Dr:4rvUI>—^U 


(a+i)‘ 


(a+i)- 


(C4) 


(^+1)  (a+i)‘ 


+/5[I>77^U  -  2I>7-i^UI> 


/ 


h6liT“( 4+i)%:«y“^ 


'B 


(C5) 


'  ^  2A+I) ^(A+I) 

(a+I ) ®( 2il)^(A+I ) ^®(a+I )  ^  ■''  A  +  I^^(A+I ) +1)  ^ ' 


°c  “  2a^+  i^(a.+i  )  ^®(a'+i  )  ®  ( ai-i )  • 


The  properties  of  the  D  and  U  operators  allow  the 
following  relationships  to  be  used  to  simplify  Eq.  (C3), 
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-  «>|0i|0>- 


^  .<'*'"l"2|0ch-“l'>2>  -  <o|oa!o>. 


^  ^  (09) 


+1  5‘=5'''l"2l°ol  5,ni>i2> 


n.+n„-l‘'^’"l"2l°6|l'''l“2>' 

(CIO) 


where 


1 1  =  •n _ i. 


°B  =  \'2^  ”2]:V;^=ryUI>- -  D - ^UI>— =5--5U 

(/!S+I)^  ^  ’  (a+I)^  (a+I)^  (A+I)^ 


(cii) 


(a+I) 


(A+I)' 


“  ^a+iT^ - - — (C12) 

0  (A+i;  (221+1)^  (a+I)  ^ 


For  the  terms  proportional  to  exp  (-^q),  the  sim¬ 
plified  expression  for  the  partition  function  to  fourth 


order  Is 


e^(/so)[i  +^[<o|D;^^|o>  -  <o|d^ud_^u|o  > 

(013) 
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+  /«V2t<0[l)^UI>^j^(0>]] . 

We  see  that  Eq.  (C13)  represents  an  expansion  to 
fourth  order,  l.e.,  through  all  terms  Involving  at  least 
two  D  and  two  U  operators,  of  the  following  exponential, 

<0|exp[-^(£^-D^^+D^U]>-^;^  (014) 


In  addition,  the  expectation  value  of  -1,^  times  the 
argument  forms  precisely  the  same  perturbation  series  for 
the  ground  state  as  was  determined  In  Chapter  III.  This 
may  be  seen  by  simply  utilizing  the  binomial  theorem  to 
obtain 


<DU>  <DIU>  ,  <DO>  , 

^  “•  O  •  •  •  f 

^  a3 


(C15) 


<0  f  DT—ryUD |0> 

(a+i)^ 

,  6<DU><DI^U> 


<W>2 


3<DUXDIU> 


•  »  *  X 
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Examination  of  Eq.  (Cl6)  shov/s  that,  to  fourth  order, 
t  represents  the  expansion  Indicated  In  Eq.  (ll4). 
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